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Preface 



Quantum Chromodynamics (QCD) is the theory thought as fundamental ^ for de- 
scribing strong interactions. The invariance under U{1) gauge transformations of 
Quantum Electro-Dynamics (QED) is extended to a non-abelian gauge one SU{3)] 
this fact implies that gluons, the intermediate bosons of strong interactions, are 
charged, and thus self-interacting. This feature makes QCD deeply different from 
the abelian theory. The study under the renormalization group of the strong cou- 
pling, hence of the j3 function, must take into account all of the gluons self-interacting 
diagrams. These terms give an opposite contribution to the screening phenomenon 
present in QED, and the consequence is a color charge growing with the distance. 
On the other hand, in the short distance regime the strong coupling goes to zero, and 
the theory displays the mechanism of asymptotic freedom (this aspect of the theory 
caused the Nobel Prize 2004 to the american physicists Gross, Wilczek and Politzer 
who had anticipated this feature in the early seventies PlISEI)- Therefore, in QCD 
we can distinguish two opposite limits: the one of high energies, where it is possible 
to apply a perturbative expansion, and the one of low energies. In the latter regime, 
the perturbative expansion becomes meaningless and in order to compute physical 
quantities we must use alternative techniques or effective models. The physical scale 
separating the two aforementioned regimes, Aqcd, is of the order of a few hundred 
MeV (from recent determinations Aqcd ~ 200 MeV). Scattering processes charac- 
terized by energies > 1 GeV have been checked to be in agreement with perturbative 
QCD predictions. In this work, we will face mainly the non-perturbative aspects of 
the theory. 

The main QCD open problem, not yet resolved in a satisfactory way, is confine- 

^In the framework of standard model. QCD can be otherwise considered as an effective theory 
deriving from more fundamental interactions (SuperSymmetry, String Theory), but in the following 
work we will not deal with such theories. 
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ment: quarks are not physical states which can be observed as quasi-particles, but 
only as bound states which correspond to the hadrons; these objects are forced to 
be color neutral. In practice, the potential energy for a couple qq grows indefinitely 
with the separation distance, so that it is impossible to separate them unless an- 
other quark-antiquark couple is created from the vacuum. Therefore, the only way to 
break a color string is to divide it into two shorter ones. Confinement deals with the 
long distances regime of the theory, and is therefore intrinsecally a non-perturbative 
phenomenon. 

Another non-perturbative feature of the theory is the dynamical breaking of 
chiral symmetry. If wc consider QCD with Nf massless fcrmions, the Lagrangian 
possesses a global symmetry SU {Nf)L(S)SU (A^/)i?. Chiral symmetry is not observed 
in nature, since quantum interactions lead to the formation of a chiral condensate 
(ff^) which spontaneously breaks the symmetry SU{Nf)L SU{Nf)R to its sub- 
group SU {Nf)v- In this case we speak of dynamical breaking since it is the dynamics 
of the interactions which is responsible for the breaking. 

Chiral symmetry breaking implies many theoretical and phenomenological con- 
sequences; moroeover, it is closely connected to the confinement problem, because 
both phenomena are related to special configurations of the gauge field. These issues 
will be discussed in the Introduction. Here we stress that confinement and chiral 
symmetry breaking are both non-perturbative effects, and we cannot deal with them 
by means of a perturbative diagrammatic expansion. We can analytically explain 
confinement only in suitable extensions of the fundamental theory. For this reason, 
a study of such behaviours starting from QCD itself is a hard goal, which has not 
been satisfactory resolved yet. On the other hand, a class of effective models are 
able to mimic the basic features of the non-perturbative regime, even if not starting 
from first principles. When dealing with effective models, which are obtained as 
simplified versions of the fundamental theory, the study with different approaches 
of a physical problem is useful to make a result stable and possibly model indepen- 
dent. In the following chapters some of these effective models will be presented; in 
particular, we will focus our attention to the ladder-QCD and Nambu-Jona-Lasinio 
(NJL) models, which have been used to carry out this specific work. 

A powerful tool to deal with the non-perturbative problem of strong interactions 
is to consider QCD on a space-temporal lattice; starting from the QCD Lagrangian, 
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and by means of numerical Monte-Carlo simulations, it should be possible to derive 
all the physical quantities related to QCD without adding new free parameters. 
The problems related to this kind of approach are mostly of numerical nature^; 
first, we have to check that finite lattice spacing does not infer the determination 
of physical observables; second, the needed power of calculus is huge, and present 
supercomputers are still unable to push the analyses to the final step. In any case, 
waiting for an improvement of resources with new generation processors, we can 
think of lattice QCD as the closest "effective model" to the fundamental theory. 

Up to now, we have considered properties of the QCD vacuum. On the other 
hand, the introduction of finite temperature and charge densities can cause the 
restoration of symmetries which arc broken in the vacuum, namely at zero tempera- 
ture and density. Actually, the restoration of chiral symmetry in the regime of high 
temperatures and/or densities is a stable feature amongst several models. Moreover, 
from the early seventies, it has been argued that the restoration of chiral symmetry 
may coincide with the transition to the deconfined phase of the theory, the QGP 
(Quark-Gluon Plasma) regime, a plasma of weakly interacting quarks and gluons. 
In this new regime, the constraint of color neutrality relaxes, and it should possible 
to think of quark and gluons as quasi-particles moving almost freely through the 
medium. 

The possibility of considering QCD beyond the vacuum state is fundamental for 
establishing a comparison with realistic experimental settings. QCD in the regime 
of high temperatures regards the physics of Relativistic Heavy-Ion collisions, where 

we can assume the formation of clusters characterized by a high amount of kinetic 
energy. The investigation of this context is stimulated by the perspective of getting 
an experimental proof of deconfinement. On the other hand, the low temperatures 
and high densities regime is relevant for studying the physics of compact stars. In 
that regime, chromomagnetic forces induce an attractive quark-quark interaction 
which can lead to an instability of the Fermi surface and to the formation of Cooper 
pairs. If we replace electrons with quarks, and phonon exchange with gluon ex- 
change, we find a complete analogy with superconductivity in ordinary materials. 
For this reason, we refer to this phenomenon as Color Superconductivity (CSC). 

^The difficulty of considering finite quark chemical potentials in lattice simulations will be 
dicussed later on in this work. 
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The physical scenario arising in these regimes is very different, and will be thor- 
oughly discussed in the Introduction. As for water we find a liquid, a solid, or a 
gaseous phase, depending on the values of the thermodynamic parameters, in QCD 
we can have a hadron phase, a QGP, or a CSC regime according to the temperature 
and densities involved: the study of the different regimes of QCD by varying the 
thermodynamic parameters defines the phase diagram of the theory. As a further 
analogy with water, we will find a critical point and a critical line of first order 
transitions. The critical point should belong to the universality class of the 3-d 
Ising model. These are spectacular proofs of universality in Physics: we find a sim- 
ilar critical behaviour for systems which differ by about ten orders of magnitude in 
temperature, and which are characterized by radically different interactions. 

In year 2000, a paper of Son and Stephanov ^| has opened a new field of research 
in literature. In the framework of a low energy chiral model, they have analyzed 
the phenomenon of pion condensation in the regime of low temperatures and high 
isospin densities. Several analyses expect kaon condensation as well in other low 
temperatures regimes. There are many reasons why this physics is intriguing. First, 
from the theoretical point of view, meson condensation can regard the physics of 
compact stars, and is closely related to color superconductivity. Second, the effect 
of a low but finite isospin chemical potential is important in heavy ion experiments 
to give a signature of the expected transition between the hadronic and the QGP 
phases. Third, following the idea of Afford, Kapustin and Wilczek jH], QCD at 
finite isospin densities can be simulated on the lattice by means of Monte-Carlo 
simulations, whereas standard techniques are unable to describe the regime of high 
baryon densities. 

These issues will be discussed more carefully in the Introduction and in Chapter 
II, but we anticipate here the motivation of the present work: the study of chiral 
symmetry restoration and meson condensation by using a microscopic model. A 
model with quarks as fundamental degrees of freedom has a wider range of applica- 
bility than low energy models. In addition, we can display chiral symmetry breaking 
and restoration as a dynamical effect. In order to carry out this project, we have 
employed the models ladder-QCD and NJL. Our goal has been to investigate a large 
slice of the space of thermodynamic parameters. In the following, a rapid outline of 
the work is presented. 
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In the Introduction, we present a general survey on the QCD phase diagram. We 
will discuss more carefully the influence of chiral symmetry breaking and restoration 
on the structure of the phase diagram, in particular in relation with the confinement 
problem. We will discuss the phases of high and low temperatures, and the related 
phenomenology. We stress the presence of a critical point in the regime of interme- 
diate temperatures and densities, since it is a promising feature to be experimentally 
highlighted. 

In Chapter II, wc will review the recent literature which has stimulated this work. 
In particular, we will be interested in QCD at finite isospin densities, where, for high 
enough values of isospin charges, the condensation of charged pions is expected. The 
formalism for describing Bose Einstein condensation in Quantum Field Theory is 
presented. In order to clarify the physical problems we have dealt with, we will 
discuss some theoretical issues concerning QCD at finite chemical potentials. 

Prom Chapter III on, the more specific research part of this work starts. Here we 
present a study of two-fiavors QCD at non vanishing isospin density which has been 
carried out in the framework of the ladder-QCD model. The steps for deriving the 
two-loops effective action, encompassing a non-trivial non-perturbative dynamics, 
are elucidated. The possibility of pion condensation at high isospin densities is 
investigated, as well as the phenomenological consequences of a non vanishing isospin 
chemical potential for heavy ions experiments. 

In Chapter IV, we will perform the same analysis of Chapter II by using the NJL 
model. Due to technical reasons, this model is simpler to work with than ladder- 
QCD, and from now on will be used as the main tool for our studies. A detailed 
investigation of the two-flavors QCD phase diagram, by varying arbitrarily the two 
light quark chemical potentials, has been done. 

In Chapter V, we will extend the analysis to the strange flavor. We have stud- 
ied the possibility of pion and kaon condensation in the space of thermodynamic 
parameters, by means of a three-flavors NJL model, by making a comparison with 
previous analyses performed within low-energy models. An original proposal for 
kaon condensation in the regime of high baryon densities has been suggested. 

In Chapter VI we will focus our attention to the critical line for chiral symme- 
try restoration. The study will be performed in a three-flavors NJL model with an 
instanton-induced interaction. Three main issues are discussed: the dependence of 



5 



the critical curve on baryon and isospin chemical potential; the sensitivity of the 
order of the chiral transition on the value of the strange quark mass; the depen- 
dence of the curve T{ijl) on the number of massless flavors Nf. Some other recent 
developments are briefly discussed too. These analyses are motivated mainly by the 
growing interest that the lattice community has shown for this subject. In this way, 
the consistency between effective models results and lattice QCD simulations can 
be checked. 

In the Conclusions, an overview of the work done is afforded. 
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Chapter 1 



Introduction: A survey of the 
QCD phase diagram 

In this introductory Chapter, we give an overall view on the properties of QCD at 
finite temperature and densities. To start with, we continue the treatment of the 
zero temperature features which will be most relevant for further applications. 

1.1 Vacuum properties of QCD 

We analyze the phenomenon of chiral symmetry breaking. If we start from the 
chiral limit, namely where the three light quarks are exactly massless, the QCD 
Lagrangian admits an invariance under chiral transformations 

U{3)l ® U{3)n = SU{3)l ® SU{3)r ® U{1)l ® U{1)r (1.1) 

Actually, in the chiral limit, left/right chirality components ^l,r = (1 ±75)^/2 
= {u, d, s) is the spinor in the Dirac and flavor space) decouple, and the global 
f/(3) flavor symmetry doubles into U{3)l (S> U{3)r. A mass term in the QCD La- 
grangian would lead to the mixing of left/right components 

m^"^ ^m{^L'^R + ^R-^L) (1.2) 

As far as the group SU{3)l SU{3)r is concerned, the symmetry reads in the 
following invariance for ^'i^i? 

*L,i? ^ exp {-iXaa^''^)-^L,R (1-3) 
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where Aa, a — 1 ... 8 are the Gell-Mann matrices (representatives of SU(3) gen- 
erators) and a^'^ the corresponding arbitrary parameters. For every multiplet of 
particles, this symmetry would imply the existence of a mirror multiplet with op- 
posite parity; this is not a phenomenon observed in nature. Actually, quantum 
interactions lead to the formation of a chiral condensate (^^) which spontaneously 
breaks the symmetry SU{3)l ® SU{3)r to its diagonal subgroup SU{3)v- The 
phenomenology which follows from the spontaneous breaking of a global continu- 
ous symmetry is characterized by the presence of a set of Goldstonc bosons, cor- 
responding to the broken generators of the symmetry group ^: in this case the 
three pions (m^r ~ 140 MeV), the four kaons {rriK ~ 500 McV) and the r] particle 
(m^ ~ 550 MeV), the meson octet, can be interpreted as the pseudo-Goldstone 
bosons corresponding to the breaking of SU{3)a ^- Their mass is little with respect 
to that of the other hadrons, and derives solely from the current mass term present 
in the QCD Lagrangian. Since the explicit breaking term is small if compared to 
the dynamical generated mass, we can still treat chiral symmetry as an approximate 
symmetry. 

The formation of a chiral condensate would break in a similar way the residual 
symmetry C/(l)i(8)C/(l)fl down to C/(l)y. This would imply the spontaneous break- 
ing of the extra U{1)a symmetry. Actually U{1)a is not a symmetry observed in 
nature, because there not exists for every hadron a partner with the same mass and 
opposite parity. Moreover a ninth particle with a mass comparable with that of the 
octet and interpretable as pseudo-Goldstone for the breaking of U{1)a is not ob- 
served. U{1)a is in fact an anomaly (it is called axial anomaly), that is a symmetry 
present at the classical level but violated by quantum interactions. This problem 
was resolved by Adlcr-Bcll-Jackiw, who noticed the singular nature of the operator 
'^{y)'if{z) when considering {y — z) — > 0. The corresponding axial current Jg*, which 

^The Goldstonc theorem applies only when a global symmetry is broken. In the case of breaking 
of a local symmetry, the degrees of freedom of the would-be massless particle give mass to the 
intermediate vector bosons. Since in particle Physics only approximate global symmetries are 
actually broken, genuine massless bosons have not been observed. Massless modes can occur 
instead in condensed matter Physics. We will encounter examples of this kind in the following of 
this work. 

^This is a commonly used, but not quite exact notation. In fact SU(i)A is not properly a group 
(the commutator of two axyal charges is a vectorial one), but the counting of Goldstone bosons is 
in agreement with this naive analysis. However, in the following of this work we will use sometimes 
this notation. 
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would be conserved at the classical level, satisfies 

d,J^ = -£^nf F^^F^^ = 2Nf Q{x) (1.4) 

where Nf is the number of fermion flavors, g the coupling constant and F^^ the 
gluon field strength tensor (F^ = e^^paF^ is the dual gluon field strength tensor); 
Qi^) ~ 32^ ^ap^afs the topological charge density. The four-divergence of the 
axial current is not vanishing but depends on the gluon configurations; this fact 
implies that this classical symmetry is explicitly broken, and that therefore no Gold- 
stone boson appears. This is the reason why the rj' {m^i ~ 1 GeV) is much heavier 
than the mesons of the octet. 

A class of gauge field solutions which give rise to a violation of axial symmetry 
are instantons. Instantons can be interpreted as classical solutions of the euclidean 
action, and are related to non trivial topological configurations (with non vanishing 
winding number) of the gauge field: actually, they connect asymptotic configurations 
characterized by different topological charges. Similar solutions must tunnel under a 
potential barrier (a double well potential is a toy description for instantons in quan- 
tum mechanics), and are characterized by a finite action; therefore, they do give 
a finite contribution to physical quantities. By reminding eq. ()1.4j) . it is clear that 
instantons break explicitly U{1)a- These solutions cannot be obtained by smooth 
perturbation of the vacuum, due to their non trivial topological structure, and are 
therefore intrinsically non-perturbative. For the case of a f/ {Nf)L®U {Nf)f> invariant 
theory, 't Hooft realized P that the coupling of quarks to such localized topological 
solutions of the pure gluonic fields produces an effective interaction between quarks 
which is still SU^Nf)^ ® SU^Nf)^ but explicitly breaks the U{1)a symmetry. This 
interaction gives rise to a mixing between different flavors; as a consequence the wave 
functions of physical isosinglet bosons 7],7]' are mixed with respect to the U{1)a in- 
variant case (in that case we would have t] ~ ^'A8\E', rj' ~ '^'Aq^E'). Simultaneously, 
the 7]' mass is pushed up by this term, in agreement with phenomenological ob- 
servations. A detailed analyses of instanton effects in QCD, and a presentation of 
the so-called instanton liquid model has been done in [Zj. However, Witten jH] and 
Veneziano [Oj argued that the solution of the U{1)a problem should not necessarily 
come from instanton solutions but rather from some other type of gluonic configu- 
rations in QCD (e.g. the ones related to confinement); this resulted in the famous 
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Witten-Veneziano formula which relates the flavor-singlet (the 77' particle) mass to 
the topological susceptibility x 



ml + ml, - 2m% ^ (1.5) 



The topological susceptibility is defined as x = ||2|e=o; ^ is the vacuum energy 
density as a function of the 9 parameter. Actually in QCD a topological term 

C0 = i9 Q{x) (1.6) 

allowed by the QCD gauge invariance, can be added to the standard QCD La- 
grangian. This parameter defines the choice of vacuum among an infinity of possible 
distinct and generally inequivalent vacua. Each 6 represents a possible true vacuum 
and there are in general an infinity of distinct theories arising from any given La- 
grangian. Due to the finite value of x ~ (170 MeV)^ in the vacuum, the r]' mass 
does not vanish even in the chiral limit (where mx = = 0). 

The spontaneous breaking of chiral symmetry produces an extra contribution 
to the effective quark mass, due to interactions. Actually, the quark condensate in 
QCD is given by the trace of the full quark propagator Sp'- 

{^^) = -i lim^^,+ TV Sf{x, y) (1.7) 

Since is a gauge invariant quantity, one can take any gauge to evaluate Sf{xi y) 
which will have a general form in the momentum space as 

Sf{p) = "^^^'l 2, (1-8) 

Due to the property Tr 7^ = 0, only a non zero B, i.e. a mass term, can induce 
a non vanishing quark condensate. According to this constituent particles descrip- 
tion, the main effect of interactions is the redefinition of an effective constituent 
mass: quarks behave as quasi-particles^ with an effective mass, given by the sum of 
the current mass (the one present in the Lagrangian) and of the correction coming 
from interactions. This fact allows us to consider the quarks inside the nucleons 
as objects characterized by a mass ~ 350 MeV, instead of by the current mass, 
which is of the order of 5 10 MeV (at the 1 GeV scale). As we can see, for the 
^Obviously, confinement prevents to directly observe quarks outside of hadrons. 
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light quarks the dynamical breaking effect is dominant with respect to the explicit 
chiral violating term. The phenomenon of mass generation cannot be seen at any 
perturbative order; actually, within perturbation theory in the chiral limit, there is 
no mixing between left handed and right handed quarks, and any diagram would 
imply a vanishing mass. Therefore, to highlight this feature it is necessary to sum 
the contribution of an infinite class of Feynman diagrams. 

For this reason, a study of such features starting from the fundamental theory re- 
quires some attention. One possibility to face this problem is to consider a non 
perturbative loop expansion of the effective action, following a method proposed 
by Cornwall- Jackiw-Tomboulis [10 , and to apply it to the QCD case. This ap- 
proach led to an effective theory, sometimes referred to in literature as ladder-QCD, 
which displays the dynamical breaking of chiral symmetry, and whose phase dia- 
gram is expected to be consistent with that of QCD. However, strong assumptions 
are necessary to write down an effective action, and in this spirit ladder-QCD must 
be considered as an effective theory, althought derived directly from first princi- 
ples. The most used microscopic model for studying chiral symmetry breaking is 
the Nambu-Jona-Lasinio model (NJL) fUEI- The model was proposed in 1960 in 
order to describe the generation of a mass gap for fermionic excitations, in analogy 
with superconductivity, and to motivate the smallness of pion mass with respect 
to baryon mass. In superconductivity, the gap is caused due to the fact that the 
attractive phonon-mediated interaction between electrons produces correlated pairs 
of electrons with opposite momenta and spin near the Fermi surface, and it takes a 
finite amount of energy to break this correlation. Excitations in a superconductor 
can be conveniently described by means of a coherent mixture of electrons and holes. 
After the proposal of QCD as the fundamental theory NJL has been used at the 
quark level. In analogy with the Fermi model for weak interactions, the exchange 
of a gluon is replaced with an effective four fermion vertex; here it is the strength 
of interaction, and not the mass of the intermediate boson, which allows a simi- 
lar identification. The attractive interaction leads to the formation of a fermionic 
condensate {^'^)', we can think of Bogoliubov fermionic quasi-particles as collective 
excitations provided with a constituent mass. The model is reliable in the regime 
where QCD interactions are "strong", at a scale ~ Aqcd, see e.g. [I3j; therefore, 
NJL can be considered as a good model for studying chiral symmetry breaking. On 
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the other hand, since gluons are not present in the model, we cannot employ it for 
describing the perturbative high-energy regime. The interaction respects the global 
symmetries of the original Lagrangian, and we can model the U{1)a breaking with 
a suitable determinantal term in flavor space ('t Hooft term). 

In this work, we will deal with the two aforementioned models. However, there 
is a larger class of models that can be employed for studying the non perturbative 
features of QCD. In the Preface we have mentioned lattice QCD as the closest model 
to the fundamental theory. Otherwise, in the spirit of the renormalization group, 
we can write down an effective theory integrating out the heavy degrees of freedom. 
We can study a theory where only the low energy degrees of freedom are taken into 
account: due to confinement, these states must necessarily be color-neutral hadrons, 
and due to the breaking of chiral symmetry, these are the pseudo-Goldstone bosons 
corresponding to the breaking of the global symmetry. The interactions of these 
low energy states are given by symmetries of QCD Lagrangian. These low-energy 
models, called chiral models, are basically non linear sigma model (the non linear 
constraint fixes the a ~ field); they are effective in the sense that they display 
chiral symmetry breaking already at the tree level. Perturbative calculation can 
be done with respect to the external momentum. Chiral models receive a large 
application in nuclear physics and in hadron spectroscopy; the lightest baryons can 
be included in the theory too. A random matrix model can be used as a QCD 
effective theory as well: the idea is to replace the matrix elements of the fermionic 
Dirac operator by Gaussian distributed random variables, while respecting all the 
global symmetries of QCD. We will often refer to these approaches in the following 
of the present work. 



1.2 QCD at high temperature and Deconfinement 

We have anticipated in the Preface the existence of a phase transition associated 
with chiral symmetry. According to the Landau classification of phase transitions, 
an order parameter may separate the low temperature-broken symmetry phase from 
the high temperature-restored symmetry phase. The fermion condensate (^^') is 
the order parameter for the chiral transition. However, as far as chiral symmetry is 
concerned, a change in the symmetry properties across the transition is associated 
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only in the chiral limit, where chiral transformations are exact invariances of the 
QCD Lagrangian. The value of the critical temperature is about 150 MeV. 

Conversely, in the pure gauge limit, namely when quarks are infinitely heavy, 
there exists a rigorous phase transition associated with deconfinement. In this case 
QCD possesses a Z{3) center symmetry {Z{N) in general for a SU (N) gauge theory), 
which is spontaneously broken by the QCD vacuum; all of the dynamics regarding 
confinement is in this case ruled by gluons ^. It is still an open question what kind 
of gluon configuration are responsible for confinement: in particular, the connection 
between confinement, chiral symmetry and U{1)a breaking must be clearly under- 
stood In the pure gauge limit QCD is strictly confining, in the sense that 
the static potential between a quark and an antiquark grows with their distance. 
For finite quark masses the discrete Z{3) symmetry is explicitly broken, and in the 
infinite distance regime the static potential becomes flat: this means the creation of 
another couple qq for breaking the color string. The order parameter related to con- 
finement is the Polyakov loop L; L is charged under the center symmetry, therefore 
L breaks Z{3). The Polyakov loop belongs to the class of Wilson loops, which 
are gauge-invariant observables, and which constitute the fundamental tools for de- 
cribing a gauge theory on the lattice. At finite temperature L{T) ~ exp{—Fqq/T) 
where Fqg is the free energy of a static qq pair for r oo. Therefore (in the m — oo 
limit) L{T) = for T < and L{T) 7^ for T > T^, where is the deconfine- 
ment temperature. Up to now, the only proof of confinement comes from numerical 
studies of the Polyakov loop on the lattice. 

In the physical case, with finite values of the quark masses, neither the chiral nor 
the Z{3) center symmetries are exact. It is however possible to define an approximate 
phase transition by studying the behaviour of the susceptibilities associated with the 
respective order parameter; for realistic values of quark masses the transition should 
be a cross-over. 

Moreover, several lattice studies have found that the peaks in the chiral and 
Polyakov loop susceptibilities are reached at the same temperature. This means 
that deconfinement and chiral symmetry restoration are strictly related to each 
other. Therefore, the increase of temperature and/of quark density beyond a line of 
critical values could cause two effects. Deconfinement implies that hadrons dissolve 

■^This is the reason why a model without gluons like NJL cannot account for confinement. 
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into their elementary constituents. Chiral symmetry restoration implies that quark 
are characterized by their current mass, rather than by the constituent one. In this 
work we will be dealing with models concerning the study of chiral symmetry, having 
in mind its connection with deconfinement. 

The experimental investigation of a QGP is a challenging enterprise. Ultra Rel- 
ativistic Heavy Ion Collisions are a tool for getting a proof of the deconfined phase: 
the energies involved allow the formation of very hot objects, with temperatures 
higher than the deconfinement one, in the state following the collision. The cooling 
of these clusters, and the consequent formation of hadronic matter in the final state, 
must however bring the memory of the deconfined phase the system has passed 
through. It is very important to identify an observable which can be identified as a 
signature of the QGP. A possible proof for the formation of a QGP is strangeness 
enhancement: in a deconfined and chirally restored plasma the ratio strange/light 
quarks is sensibly higher than in the vacuum. The latter fact reflects in the pro- 
duction of strange hadrons after hadronization. Heavy-flavored hadrons production 
and charmonium (J/ip) suppression in the final channel are probes for the QGP 
formation as well. 

In comparison with earlier expectations, in recent years the study of QCD at high 
temperature has obtained remarkable progresses. Lattice studies have provided a 
prediction for the dependence of pressure/energy density on temperature, see fig. 
11.11 For an ideal gas e/T"^ = Sp/T"^ = ess/T^ (where e is the energy density, p the 
pressure and b the Stefan-Boltzmann value) . For a weakly interacting plasma the 
discrepancy between the Stefan-Boltzmann limit must be negligible, and the system 
is almost in the ideal gas regime. Actually, perturbative methods can be applied for 
studying QCD at high temperatures and densities; goes to zero for asymptotic 
values and the system behaves like an ideal gas . In an intermediate regime where 
T,fi Aqcd the interactions can still be important. According to recent lattice 
analyses, see fig. II. H until temperature is lower than a few critical temperature Tc 
(where Tc is the critical temperature for chiral symmetry restoration), the effect of 
interaction appears to be not negligible at all, and the system seems to behave as 
a perfect liquid, with zero viscosity, instead as a perfect gas. This is the picture 
of a strongly interacting QGP (sQGP) |Tn|. In this regime, the presence of bound 
states is still allowed, as they are expected to dissolve only for T > 4 ^ 5 T^. The 
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restoration of chiral symmetry could coincide with the passage to this phase and 
not to the ideal QGP. 




100 200 300 400 500 600 

Figure 1.1: e/T^ computed from the lattice as a function of temperature. The 
discrepancy between the Stefan-Boltzmann limit is still about 20% at T = 600 MeV. 
The picture is taken from pTj . 

When referring to the experimental settings, it is important to stress that a 
thermodynamic equilibrium formalism must be adopted with a certain care. The 
collision of ultrarelativistic heavy ions causes the formation of intermediate objects, 
characterized by a huge amount of kinetic energies. The typical scale of time involved 
in these processes are of the order of a few fm/c. Since we are dealing with strong 
interactions only (the other forces can be turned off at this step of the reaction), we 
can admit that these times are large enough to favour a local equilibration; hence, 
one can introduce temperature as a useful parameter, still reminding that these 
clusters expand very rapidly. A hydrodynamic description of the expanding system 
can be given, and significant evidence of this phenomenon can be extracted from 
RHIC data [IHj. In particular, elliptic flows have been detected as stable signatures 
of this fluid behaviour ^Hj- The energy densities extracted from RHIC data seem to 
confirm the idea that a sQGP has formed (it behaves like an almost perfect fluid). 
Thermalization effects (also in the framework of a non equilibrium Quantum Field 
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Theory [201); spirit of the Boltzmann equation, can be considered to take into 

account that the system is not at the complete thermodynamic equihbrium. 

As the fluid expands, its temperature decreases: when the deconfinement tem- 
perature is reached, hadronization takes place, and hadrons are formed. Since the 
typical scale related to the hadronic regime (~ GeV) is of the same order as Aqcd, 
the hadronization mechanism is intrinsecally non-perturbative; again, some effective 
model must be used. Hadrons are the only objects detectable after the collision, 
therefore hadronic models are largely used to reproduce final molteplicities. In re- 
cent years, by following previous ideas of Fermi and Hagedorn, a simple Statistical 
Hadronization Model j^U 1221 been able to to reproduce the most part of exper- 
imental data by using a few free parameters. In this model, hadrons are considered 
as a non interacting gas, and thermodynamic parameters are fitted in order to repro- 
duce experimental data. As far as total multiplicities are concerned, hadron models 
are able to reproduce the most part of the data; hydrodynamic models are more 
appropriate to justify the geometrical distributions of final hadrons. 

1.3 The chiral transition and the critical point 

Coming back to the chiral transition, the first analyses of the phase diagram con- 
sisted in the study of the order parameter as a function of temperature and 
quark density, in the framework of a thermodynamic equilibrium formalism. To this 
end, a gran-canonical approach can be applied to a model for strong interactions. 
A group of pioneer works which contributed to the determination of the general 
structure of the phase diagram is found in refs. (231211123: those authors identified 
for the first time the presence of a tricritical point in the phase diagram, separating 
first order from continuous phase transitions along the critical line. The models 
used were NJL ^3] and ladder-QCD PHl25j. The presence of a tricritical point is a 
promising feature to highlight experimentally the passage to the deconfined phase. 
For a general treatment of the theory of tricritical points see |^. In the following, 
some general aspects of the chiral transition in QCD are reviewed. 

In QCD with two massless quarks, chiral symmetry, spontaneously broken at 
T = 0, is restored at finite temperature. The phase transition is supposed to be of 
second order j2Z| and to belong to the universality class of 0(4) spin models in 3 
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Figure 1.2: A schematic phase diagram of QCD, from i29j. The dashed hnes repre- 
sent the boundary of the phase with spontaneously broken chiral symmetry in QCD 
with two massless flavors. The point P is tricritical. The solid line with critical 
end-point E is the line of first order transitions in QCD with 2 quarks of small mass. 
The point M is the end-point of the nuclear liquid-gas transition. SC labels the 
regime of color superconductivity. 

dimensions. If the transition is second order in the chiral case, when we consider 
two light quarks with a non-zero mass there is only a smooth crossover as a function 
of T. This picture is consistent with lattice simulations and several effective models, 
however the problem is still matter of debate in the scientific community. 

On the other hand, at T=0 several models suggest (NJL [221, Ladder-QCD 
Random Matrices j2H]) that the restoration of chiral symmetry at finite /i is realized 
through a first order transition. Assuming that this is the case in real QCD, one can 
argue that there must be a tricritical point in the (T, /i) plane, where the transition 
changes smoothly from first to second order. Using the notation of [22] in fig- 11-21 
we will refer to it as P. By analogy with an ordinary critical points, where two 
distinct coexisting phases become identical, it is possible to identify a tricritical 
point where three coexisting phases become simultaneously identical. A tricritical 
point marks an end-point of three-phase coexistence. In order to highlight this in 
QCD, it is necessary to consider another dimension in the space of parameters - 
that of the quark mass m. In such a three-dimensional space, there are two surfaces 
(symmetric under m —m) of first-order phase transitions originating from the 
first-order line at m = 0. On these surfaces, or wings, with m 7^ 0, two phases 
coexist: a low density phase and a high density phase. There is no symmetry 
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distinguishing these two phases, chiral symmetry being exphcitly broken at m 7^ 0. 
Therefore, the surfaces can have a hedge which is a hne of critical points: along 
this line, the discontinuity in the order parameter between the high and low density 
phases vanishes. These wing lines start from the tricritical point. The first-order 
phase transition line can be identified as a line where three phases coexist: one 
high T and n phase and two low fi and T phases (with opposite signs of m and 
consequently of ('i'\E')). This line is called a triple line. 




Figure 1.3: Phase diagram of QCD in the space (T, /x, m) for two light flavors, from 
[2H]- In the figure only the wing with positive m is shown (the other is symmetric 
with respect to the plane m = 0). Marked lines label first-order transitions. 

The appearance of a tricritical point can be understood by considering a Landau- 
Ginzburg expansion for the order parameter of chiral symmetry breaking, = 
(cr,7r) ~ (^'\E') (a is the scalar isosinglet, tt the pseudoscalar isovector): 

f}(0) = ^ 6(02)2 + c(02)3 ^^ 9) 

being ^2(0) the appropriate effective potential. The coefficients a, b and c > are 
functions of fi and T. The second order phase transition line described by a = at 
b > becomes first order when b changes sign; therefore, the condition a = b = 
locates the tricritical point. The critical properties of this point can be derived 
from universality arguments j2Hl ISOI and the critical exponents are those of the 
three-dimensional mean field theory (apart from logarithmic correction arising from 
the renormalization group analysis [2Sj)- Actually, at finite temperature the long- 
wavelength modes live in a 3-d space, the temporal direction being compactified. 
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In |28| 123 ISni a general survey of the properties of the critical point in QCD is 
given; the great merit of those papers is to have extracted a universal picture from 
a description which is in general model-dependent. 

In real QCD with nonzero quark masses the second order phase transition be- 
comes a crossover and the tricritical point becomes a standard critical end-point of a 
first order phase transition line. This point is labeled in fig. II. 21 as E. In the expan- 
sion of eq. ()1.9|) a current mass enters as a term ~ —ma. It can be argued that the 
point E is in the universality class of the Ising model in 3 dimensions, because the a 
is the only genuine massless field at this point (the pions remain massive because of 
the explicit breaking of chiral symmetry) . It is a general property of liquid-gas type 
transitions (like that of water), not being associated with a change of the symmetry 
properties across the transition, to terminate with endpoints which belong to the 
Ising-3d universality class. Universality arguments also predict that the end-point 
E in QCD with small quark masses is shifted with respect to the tricritical point P 
towards larger fi as shown in Fig. 11.21 

The strange quark plays an important role on the location of the points P and 
E. At yU = 0, if the strange quark mass is less than a critical value rris^, the finite 
temperature transition changes from second to first order. This leads to a tricritical 
point, Ps, in the (T, rris) plane. In terms of eq. ()1.9j) the effect of decreasing rris could 
favor the coefficient b to become negative [221 ■ The instanton induced interaction for 
three massless quark would cause in the expansion a 0^ term, forcing the transition 
to be first-order even at /i = [31] . 

Obviously, the physical picture in the (T, fi) plane is that of Fig. 11.21 only for 
rris > ^s3- As rUs is reduced starting from infinity, the tricritical point P of Fig. 11.21 
moves to lower /i until, at = m^s, it reaches the T— axis and can be identified 
with the tricritical point in the (T, m^) plane The two tricritical points, P 

and Pg, are continuously connected. For rrig < rris^ the transition should remain first 
order at all nonzero /x; in this case there would not be any critical end point along 
the line for chiral symmetry restoration. In the real world, with physical values of 
quark masses, it is generally assumed that rris > ^s3, in agreement with several 
lattice results and effective model calculations. It is important to stress that the 
qualitative effect of the strange quark is to reduce the value of /ip, and thus of fiE, 
with respect to the pure two flavor case, being /xp = at rris = n^'ss- This shift 
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may be significant, since lattice studies show that the physical value of is only 
slightly higher than m^s. Therefore, a reliable study of the phase diagram must take 
the strange quark into account. In fig. 11.41 we plot a schematic picture of the phase 
diagram at = in the {mu,d, plane. 
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Figure 1.4: QCD phase diagram in the iJiu^iTng plane, from [23] ■ The finite tem- 
perature, zero chemical potential transition is here considered. For realistic values 
of quark masses, the = transition should be a crossover. 

The critical end point seems to be a rather stable feature for the QCD phase 
diagram; however its exact location in the (T, /i) plane is much more uncertain and 
model dependent. For instance, ladder-QCD provides an estimate which lies at 
higher temperatures and lower chemical potentials than that of NJL. Recent lattice 
simulations j^H EH] seem to favor a high temperature location. 

RHIC and LHC (whose realization is forecast for 2007) should be able to detect 
the critical point. Since E is a genuine critical point, one expects phenomena related 
to the enhancement of long wave-length fluctuations [22] ■ Event-by-event fluctua- 
tions in heavy ion collisions can be studied as probes of the critical point [HE] • At the 
critical point the a mass is rigorously zero even in the massive case, so there must 
be a region around E where the decay cr ^ vr^ vr^ is forbidden. This fact should 
lead to a non-thermal enhancement of low pf pion spectrum in the final channel. 
Moreover, pion properties in the final channel are strongly modified, because in this 
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region the inter-hadron interactions should be dramatically changed |37j . 



1.4 The regime of low temperatures and high den- 



After the papers [SH] and [Sn], the interest for the study of the QCD phase diagram 
has grown up considerably. It is by now a commonly accepted picture that in the 
high density and low temperature regime new phases, different from that of high 
temperatures can develop. Actually, Cooper pairs of quarks could form, namely 
pairs of fermion bound together by the chromomagnetic force (see for instance 00] )• 
The basic idea is the following: in the high density regime, where the Fermi 
momentum is higher than Aqcd, we can interpret quarks as weakly interacting quasi 
particles. Furthermore, as we know from the theory of superconductivity in ordinary 
materials, the Fermi surface is unstable in presence of an arbitrarily weak attractive 
interaction. Actually, in QCD it turns out that the effective interaction between two 
quarks mediated by the exchange of a gluon is attractive in the colour antisymmetric 
3 channel whereas it is repulsive in the symmetric 6 channel (3(8)3 = 6©3). However, 
by studying the renormalization group dependence of these effective couplings, one 
recognizes that at the Fermi surface the attraction in the 3 channel is dominant 
|41j : this may induce the formation of qq pairs. Of course, this is the dominant 
mechanism in the regime where the coupling gets weak, and multi-gluon exchange 
are suppressed. Therefore, the mechanism of coupling between quarks on the Fermi 
surface, with opposite spatial momentum, is expected. The pair behaves like a 
boson: this is the reason why we speak of color superconductivity, and of di-quark 
condensation: actually, superconductivity can be interpreted as a condensation of 
fermionic pairs. The most studied case is that of a spin zero condensate, since it 
allows a better use of the Fermi surface. Recently, the case of a spin one condensate 
has been considered In order to study the ground state, we have to consider 
the following matrix element: 



Where (p^j is a Nf x Nf matrix in flavor space (with indices i,j) and a x Nc 
matrix in color space (with indices a,b), C is the charge conjugation matrix; Dirac 
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indices are not explicitly shown. Symmetry considerations are useful to study the 
order parameter: 

• The condensate must be antisymmetric in color indices (a, b) in order to have 
attraction. 

• The condensate must be antisymmetric in Dirac indices in order to have spin 
zero. 

• Pauli principles requires antisymmetry in flavor indices, flxed the color and 
Dirac structure. 

Therefore, the structure of the ground state depends on the number of colors 
Nc, which is equal to three in the real case (even though effective theories with 
different than three are often considered), and of degenerate light flavors Nf. 

A di-quark condensate (apart from the case of two colors) is not invariant un- 
der color transformations: this implies the breakdown of color gauge invariance. 
As usual, this statement has to be interpreted with care. Local gauge invariance 
cannot really be broken, according to the Elitzur theorem Nevertheless, the 

spontaneous breaking of gauge invariance is a useful concept. We can flx the gauge, 
introduce a gauge non-invariant order parameter and study its effect on gauge in- 
variant correlation functions. The most important gauge invariant consequence of 
superconductivity is the appearance of a mass gap; gluons aquire a mass through the 
Meissner-Anderson-Higgs mechanism, and there is a gap A in fermions spectrum for 
exciting a single particle from the condensed ground state. Color superconductivity 
may also lead to the spontaneous breaking of global symmetries; broken global sym- 
metries imply the appearance of Goldstone bosons, and determine the low energy 
effective description of the system. 

The qualitative structure of these superconducting phases strongly depends on 
the densities that are at stake. For high enough densities, so that we can neglect the 
mass difference between the strange and the two light quarks (and the value itself of 
the masses), the fermionic condensates which form should be those involving u,d, s, 
as though if they were massless: this phase is said CFL {Color-Flavor Locking), 
corresponding to Nf = 3. For lower densities, the 2SC ( Two Flavor Superconductor) 
phase appears to be more favored, and only the two light flavors are involved, hence 
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Nf = 2. Other more exotic possibilities can be studied but are less likely, see e.g. 

m- 

The two situations are quite different from each other. In the case of CFL we 
have a condensate (in the simplest case where the coupling in the 3 color sector is 
considered only): 

3 

(indices L/R label left /right chirality components of the spinor) 

SU{3), SU{3)l SU{3)r0 U{1)b ® U{1)a ^ SU{3)c+l+r ^Z^^Z^ (1.12) 

We can actually believe that U{1)a is restored at high densities, therefore we 
consider an interaction which respects that symmetry. The Z2 symmetries arise 
since the condensate is invariant under a change of sign of the left- and/or right- 
handed fields. We observe that CFL implies that chiral symmetry is broken. This 
mechanism is rather unusual: the order parameter involves no coupling between 
left and right handed fermions. Nevertheless, both left and right handed flavors 
are locked to color, and because of the vectorial coupling of the gluon left handed 
components are effectively locked to right handed components. The pattern of 
symmetry breaking of CFL phase (apart from the spontaneous breaking of ^7(1)^) 
looks quite similar to that of zero density; therefore, we might expect that the low 
density phase and the high density phase could be continuously connected without a 
phase transition. Electric charge is broken but there exists a linear combination with 
a broken color generator which is still unbroken. Baryon number is broken. Since 
the color group is completely broken, all gauge particles get a mass, and all fermions 
are gapped. As far as Goldstone bosons are concerned, they are 8-1-2 corresponding 
to the broken global symmetries. 

For 2SC we have a quite different situation; the condensate is 

i^tL^'^^) = Ae^'h,, (1.13) 

Due to the antisymmetry in color the condensate must necessarily choose a direction 
in color space. The condensate breaks the color group but it does not break any 
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flavor symmetry: 

SU{3)c ® SU{2)l ® SU{2)r ® U{1)b ^ SU{2)c ® SU{2)l ® 5f/(2)R ® ?7(1) ® 

(1.14) 

therefore cliiral symmetry remains unbroken in the 2SC phase. Baryon number is 
broken, but there exists a linear combination of baryon number and of a broken 
local generator which is still unbroken, therefore no Goldstone bosons appear. Due 
to the breaking of the local symmetry SU{3)c down to SU{2)c, five of the eight 
gluons acquire mass. Up and down quarks with two out of the three colors become 
gapped. 

In order to construct an effective theory it is necessary to establish a hierarchy 
of the masses that are into play. In any case, the lowest mass degrees of freedom 
are Goldstone bosons, ungapped fermions and holes and massless gauge fields. The 
corresponding low-energy effective theory has been formulated in jlH] for CFL and 
in ITT] for 2SC phases; in |3H] a general treatment of this topic is proposed. If 
we increase energy, we must consider the degrees of freedom whose typical energy 
is of the order of the gap A. For /i ^ ^qcd the theory is perturbative, and quarks 
can be treated as almost free particles. In that regime, we can consider as light the 
particles staying in a thin shell around the Fermi sphere, and as heavy all of the 
others (in particular antiquarks for /i > 0). A high density effective theory can be 
derived integrating out all of the heavy degrees of freedom jlH]; the larger it is the 
chemical potential the better it is the approximation. 

In the case where the two fermion participating to the pair do not have equal 
momentum in magnitude, pi + p2 = 2 q 7^ 0, the condensate carries a non van- 
ishing total momentum, breaking rotational and translational invariance; this can 
be the case when we consider a system where neutrality under electric charge is 
imposed, and the two light quark chemical potentials (namely Fermi energies) are 
consequently different. The most simple form of the condensate compatible with 
a non zero q is just a plane wave, (\E'(a;)\E'(x)) ~ A e'^'^'^'^, although more general 
cristalline structures could be considered. This non homogeneous structure of the 
order parameter corresponds to the so called LOFF phase, see e.g. [SD]- If the 
difference of radii of Fermi spheres is high enough, the formation of the Cooper pair 
becomes unfavored and the condensate vanishes. 
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Recent studies have proposed the possibihty of kaon condensation in supercon- 
ducting regimes, before the pure CFL phase takes over; in that case, the non van- 
ishing value of strange quark mass plays the role of an effective chemical potential 
for strangeness |13] and favors kaon condensation. 

The study of the QCD phase diagram turns out to be therefore remarkably 
complicated by those new phases, with respect to the pioneering studies of the 
eighties. In Fig ]1.5l a qualitative picture of the situation under examination is shown. 
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Figure 1.5: Pictorial representation of QCD phase diagram, from 

Heavy ion experiments are not intended to reproduce a scenario compatible with 
superconducting effects: the most natural context this matter of studies is devoted 
to is the physics of compact stars. In astronomy, a compact star (or compact object) 
is a star that is a white dwarf, a neutron star, a quark star, or a black hole. Compact 
star is a name often used when the exact nature of the star is not known, but when 
there are suggestions that it is very massive and with a small radius, thus leaving 
the above-mentioned possibilities. For example, a typical radius of a neutron star 
can be estimated to be about 10 Km (1/100.000 the solar radius), with a mass which 
can be comparable with that of the Sun. 

When a star has exhausted all its energy (stellar death) the gas pressure of the 
hot interior cannot support anymore the weight of the star (the gravitational pull) 
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and the star collapses to a denser state: a compact star. In particular, when the star 
mass is above 1.44 solar masses (the so-called Chandrasekhar limit) enough energy 
is available from the gravitational collapse to force the combination of electrons and 
protons to form neutrons. As the stars contracts further, all the lowest neutron en- 
ergy levels are filled and the neutrons are forced into higher and higher energy levels. 
This creates an effective pressure, due to the Fermi statistics, which prevents further 
gravitational collapse, forming a neutron star. However, for masses greater than two 
or three solar masses, even neutron degeneracy can't prevent further collapse and 
it continues toward the black hole state. Quark stars are thought to represent an 
intermediate regime between the two latter, with a core of deconfined strange mat- 
ter. We can therefore believe that these astrophysical objects can be characterized 
by high enough densities to favor the appearance of superconducting phases. In any 
case, the expected temperatures in the core of such a system (~ keV) are much lower 
than the critical temperatures for melting the superconducting regimes (~ tens of 
MeV), so that the latter would not be destroyed by thermal effects PO] . 

The condensed matter regime of QCD has been indicated as responsible of inter- 
esting and up to now unresolved astrophysical problems, as gamma ray burst [SI]. 
The signals emitted from the astrophysical objects called pulsars, characterized by a 
well definite period of oscillation, can be explained as generated from a point source 
of radiation on a spinning neutron star. The occasional "glitches" observed in the 
period of oscillation can be attributed to the presence of vortices in the superfiuid 
regime within the star. Other cosmological questions, like baryon asymmetry and 
CP-violation in the universe, can be faced within this formalism. QCD in the high 
densities regime represents a bridge between historically disconnected branches of 
physics like quantum field theory, condensed matter, astrophysics and gravity. It 
is a relatively new field of application that will receive further attention as soon as 
stable signatures of compact stars will be obtained. 
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Chapter 2 

QCD at finite isospin density and 
meson condensation 

Recently (2000) Son and Stephanov have considered the possibihty of pion con- 
densation in the region of high isospin densities of the QCD phase diagram. After 
that paper, the possibihty of meson condensation, initially proposed by Migdal [^2] 5 
Kaplan and Nelson [53J and Politzer and Wise jHlI starting from the 70 's in the 
context of nuclear physics, has been reformulated in the framework of QCD effec- 
tive theories, in presence and not of superconducting effects. A flavor-asymmetrical 
composition of quark densities favours pion and/or kaon condensates. In the core 
of a neutron star, due to the charge neutrality Q = constraint, the up and down 
chemical potential are not equal, since they carry different electrical charge. The 
linear combinations fiq = (/i„ + /i(i)/2, associated with the conservation of the baryon 
number, and /i/ = — /id)/2, fixed by imposing Q = 0, are introduced. An isospin 
chemical potential fij measures the different chemical composition of u,d quarks; in 
a neutron star /i/ is negative since neutrons dominate over protons. At T = /ig = 
for > ?Tj'7r/2, according to there should be the condensed pions regime. The 
effect of a /i/ 7^ 0, besides astrophysics, could regard also the heavy ion collision 
experiments; in that context (being unlikely the phenomenon of pion condensation) 
a /i/ 7^ could induce the splitting of the light quark critical lines. Moreover, QCD 
at /i/ 7^ has a positive definite fermion determinant and can be safely simulated 
on the lattice |5^. 

In the first section of this Chapter, we consider Bose-Einstein condensation within 
the formalism of Quantum Field Theory. In the second section, we focus our at- 
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tention to QCD at finite isospin density. In the third section, we consider some 
theoretical issues concerning QCD at finite chemical potential, which are of interest 
for the present study. 

2.1 Bose-Einstein condensation in Quantum Field 
Theory 

In this section we will review the phenomenon of Bose-Einstein condensation (BEC) 
in the framework of Quantum Field Theory. Let us consider a relativistic system of 
charged, interacting bosons, described by the following Lagrangian density 

£0 = - m'r* - A($*$)' (2.1) 

where a positive A (necessary in order to have a well defined theory) describes a 
repulsive interaction. The system possesses a U{1) invariance related to the charge 
conservation. We can introduce a chemical potential fi coupled to the conserved 
current, and thus we obtain an effective Lagrangian density in Minkowski space 

C = {do + ifJ^Wido - i/i)* - (V$*)(W) - m2$*$ - A($*$)2 (2.2) 

The chemical potential can be introduced by noticing that it enters into the La- 
grangian as the zero- component of a gauge field. Therefore, we can promote the 
f/(l) global symmetry to a local one, and then replace the zero component of the 
gauge field with /i. We are allowed to consider fi > m since the theory is an inter- 
acting one for A 7^ 0. 

In order to study the ground state of the system we can study the static part of 
of ea. fjT^ 

^($) = (;x^-m2)$*$ + A($*$)2 (2.3) 

We notice that, by expanding eq. (|2.2j) and neglecting derivative terms, the chem- 
ical potential enters as a negative square mass term, thus it tends to favour the 
spontaneous breaking of the U{1) symmetry. In fact, for fi > m {fi is chosen posi- 

2 2 

tive without loss of generality) the minimum moves from ($) = to ($)^ = ^ „ 7" . 
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This state is a Bose condensate, with a charge density given by 



P 



-(/i -m ) 



(2.4) 



If we expand the action up to the second order with respect to the fields, we can 
compute the mass spectrum. We can do this by decomposing the complex field $ 

as 



and searching the energy eingenmodes as a function of X11X2 For jj, < m we have 
two modes with effective masses m ± fi. At the critical value fi = m one mode 
becomes massless; this mode represents the Goldstone Boson corresponding to the 
spontaneous breaking of the U{1) symmetry in the condensed phase. In the con- 
densed regime fi > m, there is of course a massless mode. The study of the effective 
action in the condensed phase allows to describe the system as a superfluid. 

Actually, in the language of Quantum Field Theory, the condensation of a field is 
reached when the effective mass of the field goes to zero. When the effective mass is 
zero, the energy cost for exciting a particle with zero spatial momentum is vanishing. 
Therefore, in the zero temperature limit, only the p = modes contribute to the 
ground state; all of the charge is carried by the ground state. 

Here we have considered the zero temperature, finite chemical potential regime of 
the theory; one can introduce the temperature by means of standard techniques, see 
e.g. ref. [55^. Actually, this can be done quite easily in the path integral formulation 
of the euclidean theory. It is possible to study the distribution of charge density in 
momentum space as a function of temperature. In any case when we are in the con- 
densed regime, the condensate, which is the zero momentum mode of the field, carries 
a finite fraction of the charge density of the system. As the temperature is raised 
up, some of the charge is excited out of the condensate. Eventually, the temperature 
becomes high enough to completely melt, or thermally disorder, the condensate; this 
happens at some critical temperature T^. However, in the non-relativistic limit, one 
recovers the same critical temperature as in the first-quantization formalism. The 
ground state can be interpreted as a coherent state which well reproduces a classical 
state with definite number of particles (AA^ <^ A^) in the thermodynamic limit. 

In the following of this work we will be interested in the condensation of compos- 



(2.5) 



29 



ite particles, like mesons, which cannot be described by a simple bosonic Lagrangian 
as in eq. ()2.1|) . To this end, will we deal with a class of QCD effective models. The 
condensation of pions and/or kaons is described by the non zero expectation value 
of the relative quantum field, with the consequent breaking of a phase symmetry. 

2.2 Isospin chemical potential and pion conden- 
sation 

In recent years, the study of the phase diagram of 2 flavor QCD in the space of quark 
fJ'q = ifJ'u + fi'd)/'^ = (where /ie is the baryon chemical chemical potential) and 

isospin /i/ = — /irf)/2 chemical potentials has been the subject of several papers. 
In 2000, Son and Stephanov |3] have shown that a high enough value of fij may 
cause the condensation of charged pions. By using a chiral model they showed that 
fii is basically half of the charged pions chemical potentials; since pions are zero spin 
particles we expect a phenomenon of condensation when the pion chemical potential 
is equal or greater than the pion mass (we are using the relativistic notations for 
chemical potentials). 

Chiral models are field theories built in terms of the lighter degrees of freedom 
present in the spectrum, namely the pseudo-Goldstone bosons corresponding to the 
breaking of the approximate symmetry SU{2)l ® SU{2)ji (or the larger and less 
valid SU{?>)l ® SU{3)r). The constraint of invariance under the residual SU{2)v 
of QCD produces an interaction between those light degrees of freedom, as in a 
non-linear sigma model. This interaction allows, differently from the case of free 
bosons, to have chemical potentials higher than the mass. The formation of a 
pseudoscalar condensate leads to the spontaneous breaking of parity and of a U{1) 
isospin symmetry: therefore the existence of a Goldstone mode is expected and the 
system behaves like a supefiuid. 

The chiral Lagrangian describing the low-energy dynamics is 

C = i/2Tr[9^S9^St - 2m2ReS] (2.6) 

where E = e^'^"^l^^ g SU{2) is the matrix pion field (vf are the real components of 
pion field), and /^r the pion decay constant. The theory is invariant under SU{2)v 
transformations. The introduction of an isospin chemical potential explicitly breaks 
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the symmetry down to U{l)v- In order to introduce a chemical potential a possible 
way is obtained by noticing that the chemical potential enters in the Lagrangian as 
the zero-component of a gauge field. Then one can promote the global SU (2) i ® 
SU (2)/j symmetry to be a local one, and gauge invariance fixes the way the chemical 
potential enters [S^ : 

C^jf = ^TrV^SV^S^" ^ReTrE (2.7) 

The covariant derivative is defined as 



Vo = 5„S - fijirsJ: - Srs), V,S = (2.8) 

We can see that the chemical potential breaks Lorentz invariance. If we consider 
now the static part of ()2.7|1 we can study the behaviour of S as a function of One 
can prove that, by introducing the following parametrization for S on the minimum 
of the effective potential 

S = cosa + i{Ti COS0 + T2 sin0) sina (2.9) 
the effective potential depends on a only [HZj 

Veffia) = ^V/(cos2a - 1) - f^mlcosa (2.10) 

For l/i/l lower than 771,^/2 the system lies in the ground state with a = 0, S = 1. 
When l/i/l exceeds the critical value mjr/2, the minimum of S occurs for 



cosa = 




(2.11) 



For a 7^ we have the spontaneous breaking of the symmetry associated with 
rotations in the coordinate, and consequently one Goldstone boson. We can 
interpret this phenomenon as a Bose Einstein condensation of charged pions (positive 
or negative depending on the sign of /i/). 

From 1)2.111) the chiral {uu + dd) and the pion (uj^d) condensates take the following 
values 

{uu + dd) = 2{^'^)^ac cosa (uj^d) + h.c. = 2(i|f\l/)^„c sina (2.12) 
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The chiral condensate rotates into a pion one as /i/ grows. In correspondence of 
the superfluid phase the pion field must be redefined with respect to the vacuum 
case: both the superfiuid mode and the massive excitations are suitable linear com- 
binations of the original 7r+,7r^ fields, and can be interpreted as Bogoliubov quasi- 
particles of the condensed phase. 

It is worth to stress that the ground state at > m-Tr is a Bose Einstein con- 
densate of charged pions, which are composite particles. We recall here that this 
theory is reliable as long as pions are the only particles excited by the chemical 
potential. Therefore, an upper limit for using such a theory can be identified as 
/i^r = 2/i/ ~ rrip = 770 MeV. On the other hand, in the regime of asymptotically 
high values of /i/ the description in terms of chiral models breaks down; a per- 
turbative calculation is sensible, and the result is that the ground state should be 
characterized by a non vanishing value of the field 

{dl5u) y^O (2.13) 

This condensate carries the same quantum numbers as that at smaller chemical 
potentials, but the physical interpretation is quite different. In the regime ^ 
Aqcd the condensate is a BCS condensate of d, u quarks (for positive yU/). Therefore, 
at some intermediate value of /// there should be a transition BEC-BCS, but since 
the condensate has the same quantum numbers in the two regimes the transition 
should be continuous. 

In heavy ion collision we do not expect values of /// sufficiently high to cause pion 
condensation; the attainable values of /x/ could cause a separation of the critical lines 
relative to the two light quarks, being different the respective chemical potentials 
A^n = A^q + fJ-d = fJ'q — fi'i- This aspect can have a strong phenomenological 
significance. 

There is another reason why QCD at finite isospin densities has been widely studied 
in recent years. Actually, QCD at /i/ 7^ (and fig = 0) can be simulated on the 
lattice [H], whereas the general and more physical case of /i, 7^ suffers from the 
so-called sign problem (the fermion determinant is complex in euclidean space). 
However recently some approximate solutions have been proposed to overcome this 
difficulty. When integrating out the fermionic degrees of freedom, one is left with 
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the following effective action 

Z = j Vcj) det M(0) e"^^*^^'^) (2.14) 

where are bosonic variables, and detM(0) is the fermion functional determinant, 
acting in flavor, Dirac, color and momentum-frequency spaces, stemming from a 
Lagrangian density ^fM(0)\E' after the integration over fermion fields. If detM{(j)) 
is a real-positive definite quantity, we can treat this term as a modification of the 
functional measure T>(j), and standard Monte-Carlo simulations are allowed. To 
guarantee that the measure is positive definite, generally one must have an even 
number of flavors, with det{M) real for each of them. One particular case is when 
there exists an operator P such that for each flavor 



= PMP-^ (2.15) 

In this case det M^" = det M, and therefore {det M)^ = det |Mp which is positive, 
and thus for each flavor det M is real. The full determinant is positive for an even 
number of flavors. In the case of a chemical potential real and equal for two flavors, 
relation ()2.15|) does not hold, and this brings to the well known sign problem. 

One possible way to avoid this problem is to consider an imaginary chemical 
potential instead of a real one j3]. In this case the fermion determinant is positive 
definite, but after Monte-Carlo simulations an analytical continuation is needed to 
come back to the physical case of a real /i. Some caution in this procedure is needed. 

A case where numerical simulations are possible is that of QCD at finite density 
of isospin. In that case the operator M{fi) has the following structure 

L{fi) is the Dirac operator for one flavor with chemical potential /i 
L{^y = 'j5L{—fi)'j5, and eq. ()2.15|) is satisfied by setting 

--Co) 

The relation det M{^) = \det L{^)\'^ follows straightforwardly. Therefore the 
study of QCD at finite isospin density allows to establish a direct connection between 



(2.16) 
. -L(/i) satisfies 
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analytical calculations, performed within effective models, and numerical lattice 
simulations of the basic theory; this is a reason why this argument became recently 
rather popular. 

Another case where det M is real with a chemical potential is when there exists 
an invertible operator Q such that: 

M* = QMQ-^ (2.18) 

Examples of this sort are afforded by models with four-fermion interaction, like the 
NJL model. 

2.3 General remarks on QCD at /i 7^ 

The insertion of a chemical potential breaks explicitly charge conjugation (and also 
Lorentz invariance). Futhermore, the condensation of a pseudoscalar field at /ij 7^ 
leads to the spontaneous breaking of parity; also the combined symmetry CP is 
broken. The breaking of P is not in contradiction with the Vafa-Witten theorem 
[SHI Ell, which states that, in vector-like gauge theories like QCD (i.e. theories where 
gauge invariant bare masses are possible for all fermions, which implies that both 
states of chirality appear in the theory), one cannot have the spontaneous breaking 
of parity and of vector-like global symmetries (like isospin and baryon number) . The 
theorem holds at zero temperature and chemical potentials only, and in the case of 
vanishing theta term, ^ = 0. The 9 term, see eq. (|1.6|) . is related to the so-called 
strong CP problem in QCD. A 7^ would break CP, due to instantonic effects, 
but experimental analyses put a very low upper bound on 9 {9 < 10~^). This prob- 
lem could be avoided considering a dynamical degree of freedom associated with 9 
1^ : the axion would result as a pseudo-Goldstone boson corresponding to a 
Peccei-Quinn symmetry. 

The previous arguments may have a strong cosmological significance. The conden- 
sation of a pseudoscalar field can be used as a toy model description for breaking 
CP, in the attempt of giving a unified explanation of cosmological issues by means 
of particle physics [H2] • 

The physics of QCD at finite chemical potential presents many other interesting 
features. One is related to the counting of Goldstone bosons. There is a theorem due 
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to Nielsen and Chadha jHSI stating that, in relativistically non-covariant theories, 
there is a possible mismatch between the number of broken generators and number 
of massless particles. One can distinguish two types of Goldstone bosons: those 
with the energy proportional to an even power of the momentum and those showing 
a dispersion relation which is an odd power of the momentum. The particles of the 
first type must be counted twice. QCD with a non zero chemical potential is just 
an example of a Lorentz non-invariant theory. 

This theorem can be applied for instance to the case of non zero chemical potential 
for strangeness [HI]. Let us consider a theory of interacting bosons: 

C={do + ij^Wido - - (V$*)(V$) - m^*** - X{^*^f (2.19) 

where $ is a doublet of complex fields and m is the mass term. The theory is invari- 
ant under SU{2)^U{1) transformations, and the chemical potential fi is introduced 
with respect to the U{1) charge. This simple model can mimic the behaviour of the 
four kaons, being /i a chemical potential associated with strangeness; in a realistic 
model the interaction term should be better specified. However, as far as the prob- 
lem of kaon condensation is concerned, this model encompasses much of the relevant 
features. If the chemical potential is higher than the mass m, the potential energy 
has a non trivial minimum with $ 7^ 0. This causes the breaking of SU{2) ® U{1) 
down to f/(l). A naive analysis would suggest the existence of three Goldstone 
bosons, corresponding to the three broken generators. Conversely, if one expands 
the action around the minimum of the potential, and studies the particle spectrum, 
one finds only two Goldstone bosons; this is possible because one of them has a 
quadratic dispersion relation, and the other a linear one. The counting of massless 
particles stems form the Nielsen- Chadha theorem. 

Another interesting argument concerns the behaviour of the QCD partition func- 
tion at zero temperature and finite chemical potential 

Zg(T, /ig) = j VA{det{P + m- fXg^o)f e"^^^^ (2.20) 

where Sym is the euclidean finite temperature Yang-Mills action, and we are here 
considering two degenerate fiavors with equal chemical potential fiq. The fermion 
determinant is taken for each quark fiavor, therefore it appears to the second power. 
The only place where the chemical potential enters in is the fermion determinant. 



35 



One would expect that rising up /ig from zero, at T = 0, the fermion determinant 
is modified and the system is excited from its ground state. Instead, until a critical 
value of is reached, the system still lies in the ground state, and the partition 
function remains unchanged. The critical value fig = {itln — e)/3, where tun is 
the nucleon mass and e is the nuclear binding energy, must be associated with the 
transition to the nuclear matter regime. This problem (called Silver Blaze problem 
from a Conan Doyle story 65 J appears hard to be faced by means of a functional 
integral approach. As a first step, we can consider an isospin chemical potential. 
The partition function can be written as 

Zi{T, fij) = j VA\det{p + m- mo) P e"^^^ (2.21) 

because the fermion determinant is real and positive. An explicit calculation of the 

eigenvalues of the Dirac operator at finite fij {det{p + m — /i/7o) = Y\ -^i) '^^^ 

j 

be done [20] : the result is that, even if the eigenvalues and consequently the full 
fermion determinant are modified by fij, the average over gauge configurations of 
the fermion determinant remains unchanged until /i/ < m^/2; for higher values of 
fij pion condensation occurs. Actually, the gauge configurations which contribute to 
the partition function have a functional determinant identical to that of the vacuum. 

As far as the quark chemical potential is concerned, a similar calculation shows 
that the partition function at T = is independent on fig for low chemical potentials 
inniinZI, but the determination of the critical value for fig is still lacking. 

However, the fermion determinant at /ig 7^ (and fij = 0) differs from that at 
fti ^ (and fig 7^ 0) only by a phase; the so-called Phase-Quenched approximation 
consists in taking the absolute value of the determinant, and thus removing the 
complex phase. Rigorous QCD inequalitites can be used to put some bound on 
the value of fermion determinant at finite quark chemical potential, which is the 
most interesting phenomenological situation. A 1/Nc 't Hooft expansion jni] can 
be applied to shed some light on the problem [Oni HO] • A numerical study of the 
dependence of the critical curve on fig and fij in the NJL model will be presented 
in the sixth Chapter of this work. 
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Chapter 3 

Ladder-QCD at finite isospin 
chemical potential 

In this Chapter, we present a calculation, done by using the ladder-QCD model, to 
explore the phase diagram of QCD for chiral symmetry breaking and restoration, 
and pion condensation, at finite temperature with different u, d quark chemical po- 
tentials. The steps for deriving the mean-field effective potential of the model are 
elucidated. In agreement with previous investigations, we find that a finite pion 
condensate takes place for high enough isospin chemical potential with a critical 
value at T = = which is just half of the pion mass. For /i/ smaller than 
the phase diagram in the (/ig,T) plane exhibits two first order transition lines and 
two critical ending points. The latter fact could be a relevant consequence of a 
non- vanishing /i/ in heavy ion collision experiments. We refer here to the paper [7T] . 

3.1 Introduction 

As said in the Introduction, according to different models [13 123 12H1 IHO] the QCD 
phase diagram in the plane (/ig, T) should exhibit a tricritical point. Since this point 
should be located at moderate density and temperature there is some possibility of 
observation in heavy ion experiments. An important point in heavy ion physics 
is the fact that in the experimental setting there is a non zero isospin chemical 
potential /x/. Studies at finite /x/ have been the object of several papers in the last 
few years jnOl HH UHl HH UHl HHj but mainly in the regime of low temperature and 
high baryon chemical potential. The first study of the phase diagram in the space of 
parameters (/ig,/i/,T) by employing a microscopic model was done with a Random 
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Matrix model [77]. It has been found that the first order transition hne ending at 
the tricritical point of the case fij = actually splits in two first order transition 
lines and correspondingly two crossover regions are present at low values of baryon 
chemical potential. Shortly thereafter, the existence of this splitting has also been 
shown in the context of a Nambu-Jona-Lasinio model in [THj . 

In this Chapter we will consider the effect of a finite isospin chemical potential in 
a model (ladder-QCD) where the existence of a tricritical point was shown several 
years ago [21]. The reason of doing this analysis in a model different from the 
Random Matrix or the NJL ones is due to the fact that QCD at finite baryon density 
is difficult to be studied on the lattice (however it should be noticed that recently a 
new technique has been proposed ^34^ and a first evaluation of the tricritical point 
has been given in [13]). It is therefore important to study certain features in different 
models in order to have a feeling about their universality. For instance the existence 
of a tricritical point seems to enjoy such a characteristic. Moreover, the interest 
for this topic is due to results from lattice simulations and effective theories which 
forecast the existence of a phase transition toward a regime of condensed pions at 
finite /x/ [1 [ZOl EH] ■ 

3.2 The model 

In this section we will review a model used several years ago to describe the chiral 
phase of QCD both at zero temperature [HB IH21 and at finite temperature and 
density |2 ^ l25 | l83]. This model is an approximation of QCD based on the evaluation 
of the effective potential at two-loop level and on a parametrization of the self-energy 
consistent with the OPE results. The most striking features of the model are the 
dynamical breaking of chiral symmetry and the presence of a critical point in the 
(T, /i) phase diagram, which are features generally assumed to hold in QCD. 

The effective action we evaluate is a slight modification (see for instance [23) 
of the Cornwall- Jackiw-Tomboulis action for composite operators fOl IHl] ■ We re- 
call here the major steps of this calculation. We start from the Cornwall- Jackiw- 
Tomboulis formula for the effective euclidean action F: 
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where the logarithm term is the standard one loop fermionic term. The free fermion 
inverse propagator is 

Sq^{p) = ip — m (3.2) 

being m a matrix in SU (3) flavor space; r2 is the sum of all 2PI vacuum diagrams 
with a propagator S which has to coincide with the exact fermion propagator at 
the absolute minimum of F. Thus S is the dynamical variable in this variational 
approach. However it turns out useful to trade S for 

E - -§ (3.3) 

which coincides with the fermion self-energy at the minimum of F. 

In the present model (ladder-QCD) we will make the very rough approximation 
of evaluating F2 at the lowest order. That is, we evaluate F2 at two-loops with one 
gluon exchange. The relevant Feynman diagram is given in Fig. 13.11 Of course, 
in the regime where the coupling gets small this is a good approximation, since 
other graphs are suppressed by a power of a^. On the other hand, in the low- 
energy regime this procedure could appear as rather uncertain. However, it turns 
out that this approximation works rather well phenomenologically (see for instance 
|82j): the intrinsecally non-perturbative nature of this approach is able to display the 
phenomenon of spontaneous mass generation. Therefore, at this order the effective 
action is simply given by jH2] 

F[S] = -Tr In (Sq-^ - S) - i Tr (S S) + c.t. 

= -Tr In [S^^ - E) + F2 [S] + c.t. (3.4) 

As in refs. jHH Ell, we use the following parametrization for S 

Sip) = tAip)p + B{p) + z75C(p) (3.5) 

Then, by working in the Landau gauge, it is possible to show that no renormalization 
of the wave function is required and also that the Ward identity at this order is 
satisfled by taking the free quark-gluon vertex and the free gluon propagator. Also 
the latter assumptions are in agreement with the behaviour of the theory in the UV 
regime. 
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Figure 3.1: The two-loop diagram needed to evaluate T2- 

We also consider the so-called rigid case, where the strong coupling g is consid- 
ered fixed at p"^ = M^, where M is a convenient mass scale to be determined in the 
fit procedure. In the case we would consider a behaviour of coupling, masses and 
self-energies in accordance with the renormalization group, we would automatically 
take into account the vertex perturbative corrections, i.e. ladder-diagrams, in the 
evaluation of the effective action. This the reason of the name of the model. How- 
ever, we will neglect the renormalization group corrections since they make much 
harder the evaluation of the effective potential. 

In this way the relation between the scalar and pseudoscalar contribution to S 
and the terms B, C in Eq. ()3.5|1 can be easily inverted, leading to the following 
expression for the effective action, fully expressed in terms of S (!82!) 

r[E] = r2[s] + riog[s] (3.6) 

where, separating S = + i-y^T^p 

with being the four- volume; C2 = 4/3 for A'^c = 3 is the quadratic Casimir of 
SU{?))c- Besides M, also (yf is a parameter of the model. The one-loop term is 

Tiogp] = -Tr log [S^^ - S) = -Tr log [ip - m - T.sip') - ^75Sp(p')] (3.8) 

where the scalar and pseudoscalar parts of the dynamical variable are matrices in 
SU (3) fiavor space, related to the scalar and pseudoscalar quark condensates through 
the following equation 

S,(p2) + z75Sp(p2) = (s + Z75p) (3.9) 
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The function /(p^) which contains the momentum dependence of the self-energy will 
be discussed later. The fields 



{Sab) 



(3.10) 



(Pab) 



{^ah^'^b) 



(3.11) 



will be determined by minimizing the effective action. In Eq. (|3.8p . m is the mass 
matrix in flavor space which is taken diagonal 



In the spirit of a variational approach, we must make an Ansatz on the behaviour 
of / in momentum space. Of course, this would not be necessary if we were able to 
solve exactly the theory. The function /(p^) will be chosen requiring that it goes to 
a constant for p —>■ and behaves as l/p^ (apart from logarithmic terms) for large 
euclidean values of p. The former constraint derives from the comparison with four- 
fermion interaction models, where the self energy is constant up to a cut-off scale; 
actually, those models give a reasonable description of chiral symmetry breaking in 
the IR regime. The latter constraint stems from the operator product expansion 
(OPE) results, in the deep euclidean ultraviolet limit. In practice, we want obtain 
an effective potential which traces that of the NJL model in the IR regime, plus UV 
corrections which are in agreement with the perturbative asymptotic behaviour of 
the theory. The merit of this approach is that we have started with the fundamental 
theory itself, by operating suitable approximations. 

By introducing a dimensionless variable x'^ = p^ /M"^ we will consider the follow- 
ing family of functions 



m = 




/iv(x2) _ l + x^ + + ...+x- 



(3.12) 



M l + x2 + x4 + ...+x2^-2 + a;2^ 
In the limit — > oo we get the "step" function used in jHH IH2] 



e{i-x^) + e{x^-i 



(3.13) 



M 
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The reason of choosing a test function which is the ratio of polynomial expressions 
lies in the fact that the finite temperature extension of the effective potential will 
result straightforward. 

Notice that for x ^ we get 

+ (3.14) 

Now, let us consider for simplicity the chiral limit and zero chemical potential. In 
this case it is simple to get the mass-shell condition from the one loop term in 
Eq. dSIHI) (see for instance ^) 

p' + {s)'fN = (3.15) 

where (s) is the field proportional to the scalar condensate (see Eq. ()3.10p ). If 
we want to recover, at least in the infrared regime, a free particle-like dispersion 
relation (as for instance happens in theories with four-fermion interaction), we see 
from Eqs. ()3.14|) and ()3.15p that we need N > 2. In this work we will choose N = 2. 
Notice that the choice N = 1 would lead to the following dispersion relation in the 
limit of small momenta 

p^{l-2{s)^) + M^{sf + ... = (3.16) 

This might give rise to some problems in the broken phase where the coefficient of 
can become negative. However no difficulties arise for the determination of the 

critical points where (s) ^ 0. 

We can thus evaluate explicitly the effective potential 

which is UV- finite in the chiral limit, whereas it needs to be properly renormalized 
in the massive case. We have employed the following normalization condition 

dV 



(3.18) 



In the chiral limit, this requirement is equivalent to the Adler-Dashen relation (see 
for instance j82j ) 

mlf^ = -2m (3.19) 
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where {^'^)q is the condensate for one massless flavor, and m is a current mass 
term. At leading order scales as the explicit chiral violating term m. Here we 
will require the validity of this renormalization condition at the values of the quark 
current masses. 
By defining = tria/M and 

Xa = -^(^a^a) (3.20) 

the normalization condition, using Eq. ()3.1()|1 . can be written as 

dVa 3M4 3M4 

c=--?777V (3-21) 



d{.OlaXa) min I'u'^ g^{M) 

where we have introduced the parameter 

The chemical potential and the temperature dependence are introduced following 
standard methods [HE] (see for example [23 )• particular the chemical potential is 
introduced from the very beginning via the usual substitution 'p^ {po + i^.p) (we 
are working in the euclidean space) in p appearing in the Dirac operator in Eq. ()3.8p . 
On the other hand the temperature dependence is introduced by substituting to po 
the Matsubara frequency Un = {2n + 1)itT in all the po dependent terms appearing 
in the effective action. The reason for this asymmetric treatment is that the p 
dependence in the self-energies becomes relevant only at p ^ M whereas, as we 
shall see, we will be interested in chemical potentials lower than M. We will not 
review in detail this procedure, since it is a standard technique of thermal field 
theory. However, we have to decompose the structure (p^ + (m + x f{p'^) Y) of the 
one-loop term into a product nj(P^ + ^f)-: since the latter terms can be trivially 
extended to finite temperature. Of course log(n.(p^ + zf)) = log(p^ + zf), and 
therefore the roots Zi play the role of the effective mass of some quasi-particle mode. 
The fact that some of the Zi can be complex leads to instabilities in the equation 
of state in the regime of low temperatures. This was the case corresponding to 
the choice of the = 1 parameter in the test function. In particular, for = 1 
the pressure (defined form the effective potential as p = —{V{T,^) — B), with 
B = V{T = = 0) being the bag consant) was not a monotonic, increasing 
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function of the temperature. This caused a non positive-definite energy density in 
the regime of low temperatures (e = —p + Ts + yuf^, being the entropy density 
s defined as ^^), and a corresponding negative specific heat Cy = in the 
same regime. In order to avoid these instabihties, we have chosen the test function 
corresponding to iV = 2. 



3.3 Results at finite temperature and density 

In order to obtain the effective action we need to calculate the determinant of the 
operator appearing in Eq. ()3.8p . We set to zero the strange quark chemical potential 
/is = and define /i/ = — /U(i)/2 and yUg = (/i^ + /id)/2. The operator is given by 
the following 3 ® 3 matrix in flavor space: 



^i{uJn + ifXuho + ip-l- K p/275 

-p/275 ^(^n + ipdho + ip-^ - Fd 

iuJnIo + ip-^ - Fs 

(3.23) 

where we have defined Fa = nia + f2Xa (a = u, d, s); /2 is defined in Eq. fl3.12|) . with 
N = 2. Also, = cu^ + |p p. Here p is related to the charged pion condensate 

9' ^^u^^d) - {d^,u)) (3.24) 



6M3 

We have directly set to zero the hyper-charged condensates in the strange sector, 
since we do not consider here the formation of such condensates for fig = 0. This 
point will be the subject of a further study in Chapter five. The strange sector 
thus factorizes and the determination of the phase diagram for (approximate) chi- 
ral symmetry restoration is performed by studying the behavior of the light (uu) 
and (dd) quark condensates, independently on the strange quark condensate. To 
evaluate the effective action, we perform the sum over the Matsubara frequencies 
using standard methods [85J. We have to find numerically the roots Zi which solve 
the mass-shell condition given by the vanishing of the determinant in Eq. (j3.23|) . 
Although formally straightforward, the calculation is a hard computational task. 
Actually, at each integration step on \p\, we have to solve a twentieth order alge- 
braic equation in a;„ in the u — d sector. The part relative to the strange quark is 
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obviously easier to deal with. The coefficients of this equation depend also on the 
parameters x«, Xd, P, Pu, Pd, mu, m^. 

Before discussing the results at finite density and temperature let us review how 
we fix the parameters appearing in the effective action. This is done by looking at 
T = /i = 0. The parameters that we have to fit are c, M, ruu, md, mg- These 
are obtained by using as input parameters the following physical quantities: m7r±, 
mK±, rriKO, and fx- Here we do not show the explicit expressions, but we limit 
ourselves to recall the major steps in the calculation. 

The static meson masses have been computed as second derivatives of the effective 
potential with respect to the corresponding canonical meson field. For instance for 
m^± we have 

= ^ (3.25) 

The relation between the canonical pion field vr and the field p of eq. ()3.24|l can be 
derived by studying the effective action for meson fields. The two fields are related 
to each other by a constant multiplicative factor. 

The decay constants are derived by evaluating the couplings of the canonical meson 
fields to the axial-vector currents. This can be done by considering the effective 
action as a function of space-time dependent scalar and pseudoscalar fields. Here 
we just write the expression for f^, in the chiral case, the so-called Pagels Stokar 
relation 



{2< 



/[p' + f:'Y (3-26) 



2"^ dp'^_ 

where S is the self-energy S evaluated at the minimum of the effective potential. 
A similar expression will be obtained in Chapter IV in the framework of the NJL 
model. In any case, at leading order in the current quark mass, the expressions we 
find for m^, statisfy the Adler-Dashen relation eq. ()3.19|l . 

The results of the fit are given in Table I; in Table II we show the values of the 
input experimental quantities, together with the result we get from the fit procedure. 



With these values of the parameters, we find that at T = yU = the quark 
condensate for one flavor in the chiral limit has the value 

= -(248 MeVf (3.27) 
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Parameters 


Fitted Values 


M 


529 MeV 


c 


1.0 


ruu 


4.4 MeV 


md 


6.2 MeV 


rris 


110 MeV 



Table 3.1: Fit of the parameters. 



Input parameters 


Fitted Values 


Experimental Values 




139 MeV 


139.6 MeV 


mK± 


494 MeV 


493.7 MeV 




499 MeV 


497.7 MeV 




92 MeV 


92.4 MeV 


Ik 


105 MeV 


113 MeV 



Table 3.2: Comparison between the values of the input parameters as obtained from 
the fit and the experimental results. 

whereas in the massive case 

{uu) = -(251 MeV)^ {dd) = -(253 MeV)^ (ss) = -(305 MeV)=^ (3.28) 

The values we find for the light quark condensates are in agreement with the 
commonly accepted value (^^\E') ~ —(250 MeV)^. By defining the constituent quark 
masses a la Politzer jHn] 

M^onst = W = 4ML.t) (3.29) 

we get 

Ms = 385 MeV, M^^ = 256 MeV (3.30) 

where, here and in the following, the light quarks have been taken degenerate with 
average mass rh = {niu + rnd)/2 = 5.3 MeV. For the ratio nis/rh we find a value of 
~ 21 instead of the generally assumed value ~ 25. The values we find for both the 
fitted quantities and the output parameters motivate the statement that this model 
reproduces a reasonable phenomenology. 
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In the case of T = /i = and in the chiral hmit, one also finds that in order to break 
the chiral symmetry we must have 



c< 1.37 



(3.31) 



which corresponds to 




(3.32) 



This condition is satisfied by our choice of parameters. 

When fij = it was shown that in the version of the model discussed in [21] 
there is a tricritical point in the chiral limit. The model we are presenting here is 
essentially the same as apart from some slight modifications, as for instance, 
the choice of the function f2{x). However the tricritical point is still present as it 
can be seen from Fig. 13.21 (central line), with the choice of parameters of Table I. 

A complete analysis of the full three parameter space (/ig,/i/,T) has not been 
completed yet especially in relation with the possible appearance of a hyper-charged 
condensate. A carefully study will be done in Chapter five in the framework of 
the NJL model. We have examined here the case T = fig = and we have found 
that there is a phase transition indicated by a finite pion condensate starting at 
/i/ = 70 MeV. Since in our model ^ 140 MeV we agree with the results found 
in the literature [H [771 HSl ED]- We also found that for m = the critical /i/ 
is zero; actually in that case the pion mass is exactly zero because the pion is a 
rigorous Goldstone boson. In the following we will limit our studies to the regme 
of small isospin chemical potential, say fij = 30 MeV, where we expect the pion 
condensate p to vanish. In this situation the determinant in ()3.23p factorizes and 
the effective action is given by the sum of three independent terms, one for each 
flavor. Therefore the structure of the action is analogous to the case of /i/ = 0, with 
each flavor evaluated at its own chemical potential. It follows that the two light 
flavor critical lines show the same tricritical structure exhibited from the central 
line in Fig. 13.21 for pj = 0. Consequently the phase diagram we obtain for a small, 
flxed isospin chemical potential {fij = 30 MeV) is described by the two side lines in 
Fig. 13. 2| each flavor u and d shows the same structure as the central line but with a 
sphtted chemical potential fiu,d = /^g =t f^i- Notice also that although the Figure 
extends up to zero temperature, its low temperature, high chemical potential part 
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should not be taken too seriously due to the existence of color superconductivity 
(which in the present study is not considered). 

It should also be noticed that in jHT] a study of the (/i^, T) phase diagram with 
non zero fii has been performed within a NJL model augmented by the four-fermi 
instanton interaction relevant in the case of two flavors. These authors have found 
that the coupling induced by the instanton interaction between the two flavors might 
wash completely the splitting of the flrst order transition line. This happens for 
values of the ratio of the instanton coupling to the NJL coupling of order 0.1-0.15. 
In our model all the three flavors are present, and the relevant instanton effects 
would give rise to a six-fermion contact interaction. However, the 't Hooft term has 
not be considered in the present version of the model. The six fermion term would 
give rise to a non renormalizable term, whereas the present model is expected to be 
renormalizable, being derived form QCD theory. Therefore, that term would require 
some caution to be included in the model. Due to the lack of this flavor-mixing term, 
it is still an open question whether these effects would wash out the splitting as in 
the two flavor case |87j. We also ignore the effects from color superconductivity, 
since these are present only for temperatures lower than some tens of MeV, and 
they would not regard the physics of heavy ions. What we are presenting here is 
a limited study of the phase diagram, and we restrict the analysis at small isospin 
chemical potential. A complete study of the phase diagram should take into account 
instanton effects and di-quark condensation too. However, a more detailed analysis 
will be done in the next Chapter, by using the NJL model, since it is simpler to deal 
with. 

3.4 Conclusions 

In this Chapter we have considered an approximate model of QCD (ladder-QCD) at 
flnite temperature and densities. In particular we have studied the experimentally 
important situation of a non- vanishing isospin chemical potential. This situation has 
been already explored previously by various authors and we conflrm, in particular, 
the results found in and [ZHj about the splitting of the flrst order transition line 
in the plane {fig, T), for small /i/. As pointed out in [TH] this result could be relevant 
for heavy ion physics since the flrst order transition line is split symmetrically with 
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respect to the original line at = 0. This implies a reduction of the value of the 
quark chemical potential at the tricritical point of an amount given by /ii, making 
easier the possibility of discovering it experimentally. Since it is very difficult to 
perform first principle analysis of QCD at finite baryon density, we think that it 
is important to show that certain features as the existence of the tricritical point 
and the possible sphtting of the first order transition fine are common to different 
models. This suggests that these features might possess some universal character, 
valid also for QCD. However we should stress that the result of the splitting of 
the first order transition line is strictly related to the factorization in fiavor space. 
For instance, in the two fiavor case, the four-fermi interaction due to the instanton 
effects leads to a mixing of the flavors that, if sufficiently large, might wash out the 
mixing. This point needs further analysis in the more complete scheme with three 
flavors and instanton-induced interaction. In the sixth Chapter we will perform a 
detailed study of the effect of instanton-induced interaction on the behaviour of the 
critical line, in the NJL model. We will show that, in the three fiavor chiral case, 
this coupling can drive the zero chemical potential transition from second order to 
first order; therefore the 't Hooft coupling is crucial for defining the properties of 
the chiral transition. 
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Figure 3.2: Phase diagram for chiral symmetry in the {fiq,T) plane. For /i/ = 
(central line), the cross-over transition line starts from the point Pt = (0, 135) and 
ends at the point E = (240,97). The line between E and the point P^ = (325,0) 
is the line for the first order transition with discontinuities in the (uu) and (dd) 
condensates. For fij = 30 MeV (side lines), the two cross-over transition lines start 
from the point Pt = (0, 135) and end at the points Eu = (210, 97) and Ed = (270, 97). 
The lines between Eu and the point P^ = (295, 0) and between Ed and the point 
= (355, 0) are the lines for the first order transitions with discontinuities in the 
(uu) and (dd) condensates respectively. 
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Chapter 4 



Pion condensation in a 2-fiavor 
NJL model 

We study the phases of a two-flavor Nambu- Jona-Lasinio model at finite temperature 
T, quark and isospin chemical potentials This study, related to that of the 

paper jHHj, completes the analysis performed in the previous Chapter, in the ladder- 
QCD model, where only small values of /i/ were considered. The paper jHH] was the 
first in literature to study pion condensation and the restoration of chiral symmetry 
in the scalar sector in the whole space of parameters [fig, fii,T), after the analysis 
performed within the random matrix model of jTTj . 

4.1 Introduction 

The possible formation of a pion condensate due to a finite isospin chemical potential 
/i/ has been, in recent years, the subject of several papers [U [7TJ [771 HE]- Pion 
condensation has been so far primarily investigated by means of low energy models 
based on chiral lagrangians jH IHSl EOl Ell • Although these models are well suited to 
study the phases of QCD as they have the right symmetry properties, they do not 
include the combined effect of the isospin chemical potential fij with a finite quark 
chemical potential fig in order to study the pattern of chiral symmetry breaking and 
restoration as well. In order to consider the effect of both fig and /i/ on the behaviour 
of scalar and pseudoscalar condensates, we need a model with quarks as microscopic 
degrees of freedom, and to study chiral symmetry breaking as a dynamical effect. 

The effect of a small /i/ has been investigated in ref. in the context of the 
Nambu- Jona-Lasinio model and in ladder-QCD [71]. The result is the splitting of 
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the critical curves for chiral symmetry restoration for the two hght flavors, whereas 
a full study for arbitrary /i/ has only been done in the context of a random matrix 
model [77]. Studies on the lattice have been performed at finite /i/ and /x^ = in 
refs. jZg IHOl Ea ESI El and with a finite fXg and /x/ = in refs. |M1 E3 Ei 123 ESI ■ 
In the work p7] the effect of both fig and a small /i/ has also been considered. The 
case of high /i^ and small fij has been taken into account in pHj . 

In this Chapter we extend the analysis of [ZH [ZHl where it was found that the 
first order transition line ending at the tricritical point of the case /i/ = actually 
splits into two first order transition lines and correspondingly two crossover regions 
are present at low values of baryon chemical potential. In particular we will be 
working in the context of a NJL model with a form factor included such as to imply 
a UV decrease of the fermion self-energy compatible with the operator product 
expansion. In the present version of the model, we do not consider instanton effects. 
As mentioned in the previous Chapter, it should also be noticed that the coupling 
induced by the instanton interaction between the light flavors might completely 
wash out the splitting of the first order transition line. This effect was studied in 
ref. j82| in a NJL model. To shed some light on this controversial point it should be 
also important to consider the explicit SU{2)v breaking due to the mass asymmetry 
between u and d quarks, and the possible restoration of U{1)a symmetry, due to the 
damping of instanton effects in the regime of high temperatures and/or densities. 
These aspects will be discussed, although partially, in the sixth Chapter of this work. 

In section II we summarize the relevant features of the NJL model we have 
here considered, with an isospin charge included. The expression for the one-loop 
effective potential and the values of the fit parameters are reported. In section III 
we describe the bosonization procedure, which is necessary to express the effective 
action in terms of the composite fields. In section IV we discuss the several phases 
of the model, together with the corresponding symmetries, by studying the behavior 
of the scalar and pion condensates with respect to the thermodynamic parameters 
T, fig, fij (or flu, fid)- Results are shown for growing temperatures, starting from zero 
up to temperatures above that of the critical ending point. In section V we present 
a calculation which confirms the analytical agreement between fif and half of the 
pion mass, in the framework of the NJL model. Section VI is devoted to conclusions. 
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4.2 The model 



Our purpose is to explore the structure of the phase diagram for chiral symmetry 
breaking and pion condensation in QCD at finite temperature and quark densities, 
by using a microscopic model with quark degrees of freedom. This task has been 
accomplished, up to now, in the context of a random matrix model simulating QCD 
with two flavors [77j and, in the case of small differences between the u and d quark 
chemical potentials, also in the Nambu-Jona-Lasinio model (NJL) [ZH! and in ladder- 
QCD [71^. 

One motivation for using a model with quarks as microscopic degrees of freedom 
is that it gives us the possibility of studying chiral symmetry breaking and pion 
condensation at finite isospin and baryon chemical potentials; this study is not 
possible within effective chiral models. 

In ref. [ZH|, the authors made use of the NJL model with a suitable form factor 
included in the quark self-energy to mimic asymptotic freedom pSl l3Uj. This version 
of the NJL model turns out to be very close to ladder-QCD as developed in refs. [23 
|H2] where the momentum dependence of the quark self-energy is consistently dictated 
by the study of the Schwinger-Dyson equation within a variational approach (see 
the previous references for details). However, although ladder-QCD is a covariant 
and self-consistent approach, the dependence on the four-momentum of the quark 
self-energy makes the numerical computation of the one-loop effective potential with 
finite quark densities much more onerous than in the NJL case, where the quark 
self-energies depend only on the three-momentum. For this reason, in the present 
work we study the NJL model. It is reasonable to expect that when employing 
ladder-QCD, the resulting physical picture does not considerably differ from that of 
the NJL model. This has been the case in previous applications too jT^ 123113 • 

As already said in the Introduction we are not going to consider the effects of 
di-fermion condensation. Therefore our results can be considered valid only outside 
the region of the color superconductive phase, which roughly should take place for 
fig > 400 ~ 500 MeV and T < 50 MeV. For the same reason we will not consider 
values of fiu or /i^ higher than 400 ~ 500 MeV, where other di-fermion condensates 
might arise (see for instance ref. |100j ). 

Let us now consider the Lagrangian of the NJL model, in Minkowski space, for 
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two flavors u, d with the same mass m but with different chemical potentials and 

•£ -^0 ~l~ ■£/! ^" ^int 

3 



^id^ - m^'^ + '^^ + —^^U^Ta^Y+[^i'^->Ta^Y (4.1 

a=0 



where \&=(^], A=(^" lis the matrix of chemical potentials and Ta 



d J' \ fid 

a = 0, 1, 2, 3, is the set of the three Pauli matrices plus the identity. The interaction 

term of eq. ()4.ip is invariant under SU{2)i ® SU(2)r ® f/(l)y, since we want to 

mimic the spontaneous breaking of chiral symmetry. In eq. ()4.1|) we recover all the 

SU{2)a invariant combinations, for instance 

(i|r^)2 + (i|r^^5f^)2 (4.2) 

Actually the a particle (scalar with isospin / = 0) is the chiral partner of the three 
pions (pseudoscalar with / = 1). 

The interaction term of eq. ()4.H) is also invariant under U{1)a transformations since 
we are not considering here the instanton-induced interaction which would explicitly 
break U{1)a- Therefore, in eq. (j4.H) we find terms like 

+ (^^Z75^)' (4.3) 

which are invariant under U{1)a- Actually, if U{1)a symmetry was observed in 
nature, there would be a neutral, 1=0 pseudoscalar meson with the same mass as 
the cr. Phenomenologically, this turns not to be the case. 

We note that we can express either by using the variables fiu, fJ'd or the two 
combinations fig = fis/^ = ^" ^ ^'^ and fij = ^ which couple to one third 
of the baryon charge density and to the third component of isospin respectively 

= fig + fii ^Vs^ (4.4) 

To study whether a pion condensate shows up, we need to evaluate the effective 
potential. This is obtained by using the standard technique to introduce bosonic 
(collective) variables through the Hubbard- Stratonovich transformation and by in- 
tegrating out the fermion fields in the generating functional. This procedure will be 
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discussed in the next section. However, the effective potential that we have consid- 
ered is not directly obtained from the Lagrangian in Eq. ()4.1|) . To mimic asymptotic 
freedom we have to include "by hand" a form factor F{p) coupled to each fermionic 
leg as in ref . |HH] , and thus we follow the same procedure as in refs. jHOl EH] • We 
could derive this momentum dependence of the fermionic fields by starting from a 
non-local action 30j. Each condensate involving a quark and an anti-quark field 
couples to -F(p)^. The result is a one-loop effective potential which generalizes that 
of the theory described by the Lagrangian in Eq. (j4.H) . and which reduces to it in 
the limit of a constant form factor F{p) = 1 (in that cut-off scale on UV 

momenta has to be introduced): 

V = ^ixl + Xl + V) + V^^g (4.5) 
%g- ^ A p 75 h, J 



(4.6) 



hf = {iujn + Pfho - p-7 - (m + F^(p) A x/) 



where Un are the Matsubara frequencies and the dimensionless fields Xf aiid p are 
connected to the condensates by the following relations 



p = - G 



A 

(4.7) 

A 



and are variationally determined at the absolute minimum of the effective potential. 

In the previous equations, A is a mass scale appearing in the form factor -F(p^) = 
A2 



A2 + p2 

It is worth noticing that the one-loop effective potential in Eq. ()4.fij) has the same 
expression as the one derived in the previous Chapter within ladder-QCD, where 
the form factor, in place of F'^{p) in Eq. ()4.6|) . was a function introduced after 
the study of the one-loop Schwinger- Dyson equation for the quark self-energy. The 
only difference is that depends on the three-momentum whereas the quoted test 
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function depends on the four-momentum; in addition the two asymptotic behaviors 
are different (~ in the test function of ref. [7T] and ~ 1/p ^ in -F^ of Eq. 
()4.6|) ). Otherwise the two effective potentials would be identical. This observation 
also explains the reason why we have adopted the NJL model instead of ladder-QCD 
to generalize the analysis of ref. jTTj at high isospin chemical potentials: in the NJL 
model the numerical analysis is much simpler. 

To fix the free parameters of the model, which are A, the average current quarks 
mass m = {rriu + md)/2 and the coupling G = g/A"^, we work at zero temperature 
and quark densities. We first choose the mass scale A within the range A ~ 500 -7- 
600 MeV. Then we determine the strength of the coupling g and the mass parameter 
m by requiring a light quark condensate of the order (^j\E'/) ~ —(200 MeV)^ and 
a static pion mass m^r — 140 MeV. The latter quantity is evaluated, as in Chapter 
three, through the curvature of the effective potential in the direction of the canonical 
pion field [H2] (having fixed /,r at its experimental value). Actually, the relation 
between canonical pion field and p involves a decay constant 

The output parameters are the following 

A = 580 MeV; g = 7; m = 4.5 MeV (4.8) 

With these values we obtain a condensate (^^/\E'/) = —(172 MeV)^ and a constituent 
quark mass Mf = 428 MeV (defined as in jTSjEOj). The critical isospin chemical 
potential at zero temperature turns out to be pf = 89 MeV. The discrepancy of 
about 25% (we recall that the expected value of fif would be 171^^/2 ~ 70 MeV) 
is entirely due to the approximate fit procedure; the fact is that the inclusion of 
a non-covariant form factor -F(p ^) complicates remarkably the calculation of f-^. 
Actually, in our previous work based on ladder-QCD [ZI], where f-„ was consistently 
calculated within the model, the result we obtained was pf = m.,^/2. However, in 
section 1331 we will show that the equivalence between m.^12 and /i^ is an analytical 
result in the framework of the "standard" NJL model, i.e. without the form factor. 

4.3 The procedure of Bosonization 

In this section we review the standard procedure of bosonization, which is useful for 
introducing bosonic degrees of freedom in the effective potential, by integrating out 
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the fermionic variables. It is particularly important for describing the condensation 
of fermionic couples, which behave as a composite bosonic field. We will refer here 
to the SU{2) flavor case. The generalization to SU{3) is straightforward, and will 
not be explicitly described in the next Chapter. We consider here the case of the 
standard NJL model, i.e. without form factor and with a cut-off A on spatial 
momentum. 

We start from the Lagrangian density eq. ()4.H) . and we write the generating 
functional for the Green's functions in the path-integral formulation: 



Z[r],f]] = j T){'^,'^)exp{i j d^x - m)^ + 



G 



a=0 



} 



A' / V{^, ^)V{aa, 7ra)exp{i / d'^x 



(4.9) 



} (4.10) 



where 7],7] are fermionic sources coupled to we can pass from eq. ()4.1()|l to 

eq. ()4.9j) by performing the gaussian functional integral over cXa, Ha fields. 

We have started from a four-fermion interaction and we obtain a Yukawa in- 
teraction with auxiliary fields Ua, iTa- If we study the equations of motions for 



the auxiliary fields deriving from the action in eq. ()4.10jl . we obtain the following 
relations: 



aa = -G^Ta^ (4.11) 
TTa = -G'^r„Z75^ (4.12) 

When dealing with eq. ()4.1()j] . it is possible to integrate over fermionic fields; by 
setting to zero the fermionic sources, and by introducing bosonic sources J^, Jn 
coupled to auxiliary fields, we obtain a generating functional expressed in terms of 
cra,na fields: 
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■ Det 



(4.13) 



id-m + 7oA - Ta<Ta - Tai'^^'Ka 

W j Viaa, na)exp{t j d^x ^^((t^ + vr^) + + J>, 

exp{ / ? tr In - m + 70A - TaCTa - rai757ra }} (4.14) 



the functional determinant Det in eq. ()4.13|l arises after the functional integration 
over fermions, and we have used the identity In Det = Tr In to pass from eq. (j4.13p 
to eg. ()4.14j) . The operators Det and Tr act on Dirac, flavor, colour and space-time 
spaces; tr acts on Dirac, flavor and colour spaces only. 

For vanishing sources we obtain the effective action for auxiliary fields in the ap- 
proximation of one fermionic loop: 



^{crl + T^l)-i tr In \id - m + 70A - (a, + i75vr,)r„] } (4.15) 



When we are interested in the properties of the ground state only, in particular 
in the dynamical breaking of the chiral symmetry, we can consider constant fields 
<ya{x) = cFai T^a{x) = tTq. Actually the minima of the effective action correspond to 
constant fields, in order to ensure the invariance of the vacuum under rotations and 
translations. In this case we get the one loop effective potential: 



V 



^ 5Z(^a + VT^) + i tr In id-m + -^qA - (cr„ + i'y57ia)Ta 



(4.16) 



In order to obtain the final expression ()4.5p of the effective potential, it is necessary to 
replace cTq, Tr^ fields with uu, dd, u'y^d. This can be done by expanding the relations 
()4.1Hl . ()4.12j) in flavor space. Charged scalar (~ ^ri\E') and neutral pseudoscalar 
(~ i|'i75\E') combinations can be set to zero from the beginning, whenever they 
are not expected to get a non-vanishing expectation value. The expectation values 
(uu), (dd), (u'y^d) are given by the location of minima of the effective potential in 
the space of corresponding fields. 
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The approximation of taking constant fields is valid at the mean-field level; we 
consider the propagation of constituent quarks as quasi-particles with an effective 
mass 

Mj = m - 2G{^f'^f) (4.17) 

where / is a flavor index. In this model without flavor mixing, for each flavor 
the contribution to the constituent mass arising from interactions is proportional 
to the corresponding scalar condensate. In the chiral case (or for current masses 
~ 10 MeV), a constituent mass ~ 350 MeV is provided by means of interactions 
only. This is the phenomenon of mass generation. 

In the following of this work we will perform calculations at the mean-field level, and 
therefore we study the effective potential as the fundamental tool for defining the 
vacuum properties. However, we could go beyond mean-field and consider space- 
time dependent fields 



(Ta{x) = {aa) + Sa{x) (4.18) 
7^a{x) = {7Ca)+Paix) (4.19) 

The inclusion of fluctuations s, p is important for defining the properties of bound 
states, both in the vacuum and at finite temperature and densities. For instance, 
if we want to compute the dynamical mass of pions we have just to expand the 
effective action up to the second order for the fields 7!'a{x), then operate the functional 
derivative 

, , //^ , = G,\x - y) (4.20) 

where Gab{,x) is the boson propagator; obviously G.„{p)^^ = defines the on-shell 
condition. By doing this, we can prove that an isospin chemical potential enters in 
the effective pion mass just like as the half of the charged pion chemical potential, 
see e.g. |101j : 

■m^+ = uIt, — 2/i/ ; m^r- = m.,, + 2/i/ (4-21) 

The latter expressions hold when we are in the phase with zero pion condensates, 
namely for < m^/2. Actually, if we go in momentum space and consider 
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momentum-dependent fields, we are taking into account their kinetic terms. There- 
fore, we are considering the fields a a, 7!' a as dynamical, propagating degrees of free- 
dom, which can give a sensible contribution to the thermodynamic functions. How- 
ever, as for the calculation of the effective action as a function of the fluctuations 
can be done by considering the NJL model without the form factor. 

4.4 Phase diagram for chiral symmetry breaking 
and pion condensation 

In order to discuss the structure of the phase diagram, it is worth summarizing 
the symmetries of the Lagrangian density in Eq. (jUT)). Both Cq and Cint are 
SUl{2) ^SUji{2) ® t/y (1) Cg)f/J^(l) invariant. The symmetry is reduced by the mass 
term C^n to SUy{2) (g) f/y(l) and further reduced from the term which selects 
a direction in the isospin space, as is evident from Eq. ()4.4|1 . unless /i„ = fid and 
thus fij = 0. The remaining symmetry can be expressed either as Uy{l) ® Uy{l) or 
Uy{1) <S) Uy{l), depending on the basis of the fields that we are choosing. 
The baryon number symmetry Uy{l) is dynamically respected, whereas a non van- 
ishing v.e.v. of the p field defined in Eq. fl4.7|) may appear, which dynamically 
breaks Uy{l). This implies the appearance of a Goldstone mode, which is either the 
charged tt^ or tt~ at the threshold, depending on the sign of fij, whereas the other 
two pions are massive modes. The nice feature is that this is a genuine Goldstone 
boson, rigorously massless even in presence of a finite current mass. 
As far as the scalar condensates Xu,Xd are concerned (see Eq. ()4.7|l ). they do not 
break any symmetry when a non zero current mass is considered. However, since 
the mass term is small, their value is almost entirely due to the approximate sponta- 
neous breaking of chiral symmetry. Consequently we distinguish regions where the 
dynamical effect is relevant, from regions where the scalar condensates are of order 
~ m/A, namely where only the effect of the explicit breaking of chiral symmetry 
survives. 

The study of the various phases has been performed numerically by minimizing the 
one-loop effective potential. We start by showing the results in the plane, 
for fixed values of the temperature. Different regions are labelled, as in ref. \77\ . 
by the symbol of the field which acquires a non vanishing v.e.v. due to dynamical 
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effects, wfiereas tlie otlier fields are vanisliing (p), or of tlie order ^ m/A (%„ and/or 
Xd)- 

Solid lines refer to discontinuous transitions and dashed lines to continuous ones. 
However, we recall that strictly speaking only the lines surrounding regions with a 
non vanishing field p refer to genuine phase transitions, associated with the breaking 
and restoration of the Uy{l) symmetry. 

4.4.1 Zero temperature 

In Fig. 14.11 we show the phase diagram in the (pu,^^) plane at zero temperature. 
Let us start from the vacuum at T = /i„ = /i^ = 0, at the center of the picture. 
Here, in the chiral limit, the pions are the Goldstone bosons associated with the 
spontaneous breaking of SU{2)l ® SU{2)fj. We have chosen these variables in or- 
der to compare the structure of the phase diagram with that obtained in ref. [77] . 
However, if we want to recover known results given in terms of the quark chemical 
potential, we have to move from the center along the diagonal at fiu = Pd and thus 
at /i/ = 0, and increase the absolute value of Pq. At ^'^l = |^^| = 293 MeV we 

meet the approximate restoration of chiral symmetry due to the sudden jump of the 
condensates of the two (degenerate) quarks to values of order ~ m/A, which is a 
discontinuous transition. The same thing happens by moving along lines parallel to 
the main diagonal in the region labelled by Xu, Xdi enclosed between the two dashed 
lines at ^ ^'^^ = = 89 MeV (see also Fig. 14.21 where it is shown that the 
critical value of /i/ at T = is independent on pq) and by varying \fiq\. The regions 
in the top-right and bottom-left corners of Fig. 14.11 thus have the Uy{l) (8> Uy{l) 
symmetry of C with p = and Xu,Xd of order ~ m/A. 

By moving from the center along the diagonal at p„ = — p^ (and thus Pg = 0) or 
parallel to it, when crossing one of the two dashed lines at |p/| =89 MeV (we have 
already discussed the origin of the discrepancy between this value and half of the 
pion mass in the model), the absolute minimum of the effective potential starts to 
rotate along the p direction. We thus have a continuous breaking of Uyil) and a 
second order phase transition with one Goldstone mode which is at the threshold, 
right along the dashed line, either the 7r+ (in the upper part of the diagram) or the 
7c~ (in the lower part). In the chiral limit these two dashed lines merge together in 
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coincidence with the diagonal at /i/ = as the pion becomes massless in this hmit 
and the rotation is sudden, giving first order phase transitions for pion condensa- 
tion. In this case there are two Goldstone bosons associated with the spontaneous 
breaking of two U{1) symmetry groups {U^{1) ®Uy{l)) [Sn]- We recall here that in 
the present application we are not considering the explicit U{1)a violation; there- 
fore, in a more realistic NJL model, with 't Hooft determinant and U{1)a breaking, 
only one Goldstone boson, corresponding to the spontaneous f/y(l) breaking, would 
exist. 

Coming back to the massive case, and still with reference to Fig. 14.11 we conclude 
that by considering not too large and by growing we find a second order 
phase transition with the rotation of the scalar condensates into the pseudoscalar, 
namely we are faced with pion condensation in a relatively simple picture. A dif- 
ference with ref. jTTj is that we do not find the vanishing of p for values of \fii\ 
high with respect to the pion mass, but still sufficiently low to avoid considering 
superconductive phases (actually, for very low \pg\ this transition would occur in 
the present model for ~ 1 GeV). 

To explore the possibility of multiple phase transitions and thus of a richer phe- 
nomenology, we need to grow |/iq| as for instance we do in the case described in Fig. 
14.31 where we follow the path of the solid line a in Fig. 14.11 at fig = 170 MeV for 
growing fij > 0. The fields Xu and Xd are almost degenerate, both in the region 
of the approximate dynamical breaking of chiral symmetry (below fij = 89 MeV) 
and in the region of spontaneous breaking of Uy{l), where they rotate into the p 
field. Then, when the line a in Fig. 14. II crosses the solid line surrounding the region 
labelled by Xd, we see that p suddenly jumps to zero with the restoration of the 
Uy{l) group and that Xu and Xd split. Actually the latter suddenly acquires a value 
due to the dynamical breaking of chiral symmetry whereas Xu undergoes a further 
decrease and remains of order ~ m/A. 

In Fig. 14.41 we plot the behavior of the scalar condensates Xu and Xd vs. pi at 
Pq = 210 MeV, namely by following the path described by the solid line b in Fig. 
14.11 We see that we never cross the region with p ^ and that we simply pass from 
a region where the dynamical effect of the breaking of chiral symmetry is entirely 
due to a large value of Xu and Xd of order ~ m/A (at large negative pi and small 
Pu), to a region where this effect manifests itself with a large value of Xd and Xu of 
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order m/A (at large positive fij and small fid)- The region in between has almost 
degenerate and both large Xu and Xd- Pion condensation does not occur for this 
value of /ig (see also Fig. I4.2|) . 

Finally, in Fig. 14.51 we plot the behavior of the condensates at fixed fiu = 200 MeV 
vs. fid (see again Fig. 14.11) . The rotation of the pion condensate into the scalar 
ones occurs when the vertical line at fiu = 200 MeV meets the dashed line at 
fij = 89 MeV, which happens for fi^ of few MeV. Then, when fi^ has sufficiently 
increased, Xd falls to a small value of order ~ m/A with a discontinuous transition, 
whereas Xu remains constant at its large value. 

4.4.2 Finite temperature 

The evolution of the phase diagram for growing temperatures is easily understood 
as far as the regions with p = are concerned. Actually, in this case the effective 
potential at the minimum is the sum of two independent terms, one for each flavor, 
and the results are straightforwardly given through the analysis of chiral symmetry 
breaking and restoration for a single flavor at finite temperature and chemical po- 
tential (see for instance refs. [23IHni)- In Fig. 14.61 we show the phase diagram at 
zero, or small isospin chemical potential (see also refs. [ZDIZHI)- From this picture 
we see that moving along any of the critical lines of chiral symmetry restoration at 
fixed /i/, the critical value of the quark chemical potential Pg decreases for growing 
temperatures. Furthermore, for temperatures below that of the critical ending point 
E, T < T{E) = 85 MeV, the transitions are always discontinuous whereas they be- 
come cross-over transitions for T > T{E). Consequently the regions labelled by Xu 
and/or Xd in Fig- 14. II shrink when growing T and their rectilinear sides become lines 
of cross-over transitions for T > T{E) (see Fig. 14.11 and Figs. 14.71 14.81 where we 
plot the phase diagrams in the {pu^ Pd) plane at T = 60 MeV and T = 140 MeV, 
which is respectively below and above T{E)). The new feature concerns the regions 
with p 7^ 0, which also reduce their size for growing T, whereas the order of the 
transitions starts to change from first to second, beginning from the critical points 
at highest \pi\, until they reach the points of the boundaries which coincide with 
those of the regions labelled by Xu or Xd (see Fig. I4.8|l . Also the length of the 
curves of second order phase transitions to pion condensation at fixed values of \fiq\ 
sensibly reduces from low temperatures to high temperatures (see again Fig. 14.11 
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and Figs. 14.71 14. 8|) . A similar behavior, for high T, is found in ref. [77]. For 

T > T{Pt) = 174 MeV, which is the cross-over temperature at zero chemical po- 
tential (see Fig. 14. 6p . all these regions disappear from the phase diagram, which is 
thus characterized by p = and Xu, Xd ~ m/A. 

The behavior of the scalar and pseudoscalar condensates at T = 60 MeV are much 
similar to those at T = 0. As an example we plot, in Fig. 14.91 the condensates at 
fiu = 200 MeV and T = 60 MeV vs. fid (compare with Fig. 14. 5|) . The situation 
is different if we consider temperatures above T{E) = 85 MeV. For instance, at 
T = 140 MeV, we see from Fig. 14. 8| that the structure of the phase diagram is only 
slightly modified with respect to the case of two independent flavors which undergo 
cross-over phase transitions at sufficiently high values of their own chemical poten- 
tials (actually the region of pion condensation has sensibly reduced with respect to 
Figs. 14.11 14. 7|) . The phase transition associated with the spontaneous breaking of 
f/y(l) can be both second or first order, depending on the path followed. In Fig. 
14. 8( the solid line a refers to a path at /ig = vs. /i/ > where the transition to pion 
condensation is continuous. The behavior of the condensates relative to this path is 
plotted in Fig. 14.101 In Fig. 14.111 we plot the scalar and pseudoscalar condensates 
at = 200 MeV and T = 140 MeV vs. /i^. 

In Fig. 14.121 we plot the value of the critical isospin chemical potential fif beyond 
which a pion condensate forms vs. temperature T at zero quark chemical potential 
{iq. The growth of is easily understood since the pion mass (which should be 
twice fif) is expected to grow near the critical temperature for chiral symmetry 
restoration |l()2j : actually, in the chirally-restored phase pion cannot be interpreted 
anymore as a pseudo-Goldstone boson. On the other hand, no phase transition 
to pion condensation is found above the cross-over temperature for chiral symme- 
try restoration. Thus the line of critical values ends at T = 174 MeV. This is a 
somehow empirical result. There must be a deeper connection between the critical 
temperature for an interacting Bose-Einstein condensate of pions and the critical 
temperature for chiral symmetry restoration. 

The pion condensate is also expected to decrease for growing temperatures; in Fig. 
14.1H1 we show p vs. T at /ig = and fij = 200 MeV. Similar behaviors are obtained 
for fixed values of /ig and /x/. 
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4.5 Analytical equivalence between /i^ and 

In this section we present a calculation which shows the analytical agreement be- 
tween /i^ and at /i^ = T = 0, by making a (1/m) expansion in the NJL 
model. We will consider the case of two light quarks with the same mass m. Let 
us first consider the T = /x = effective potential for one flavor, with vanishing 
pseudoscalar condensate 



A 2 o roo 

^ = - - / dp pWp' + im + F^Ax)' (4.22) 

8G TT^ Jq 

Here x is related to the scalar condensate through eg. ()4.7j) . and F"^ = F^p^y. In 
the chiral case we have 

A2 o roc I 

^=8G^'~^l dpp'^p' + F^A^xl (4.23) 

Xo standing for the dimensionless field in the chiral case. We find the gap equation 
by imposing the stationarity condition dV/dx L =0 

Tdp P'^^^^— = (4.24) 

Here we are interested in minima of the effective potential characterized by Xo 7^ 0, 
namely in solutions of 

13/"°° 

4;-47 / dpp^ , = (4.25) 

In the massive case, the gap equation becomes 

Ay 3 /"^ , 5 (m + AFMF^ 
AG Jo "^"^ + (m + KF'^xY 

If we consider low enough values of masses, we can separate the scalar field into the 
chiral part plus a mass correction 



X ~ Xo + Xi AXi = a m (4.27) 
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a plays the role of cliiral susceptibility a = d{Ax)/dm\m=o-- After some tedious but 
trivial calculations which follow from eqs. ()4.25p . ()4.26|) we find for a 



a 



7r2 , 



FY 



dp 



(p2 + A2F4^2)3/2 



/ 



1 

AG 



p^F^ 



dp 



(p2 + ^2^4^2)3/2 



(4.28) 



We will use this result later on. Now, we can study the effective potential for two 
degenerate flavors with mass m, with a non-vanishing pion condensate, ai fiq = T = 
and fii ^ (we are therefore considering the regime /i„ = —fid)- In this particular 
case the expression for the effective potential turns out to be quite compact: 



V 



4G^^ 



p') - — / dp p" .J{E + /xj)2 + F^AV^ + ^{E - ixiY + F4A2p' 

(4.29) 

where Xu = Xd = X (th^ two fields are equal for opposite values of the chemical 
potential), E = ^Jp^ + (m + KF'^x)'^- Now we should find the minimum of the 
effective potential with respect to the variables Xi P'l order to do this, we will 
make the reasonable assumption of taking the scalar condensate fixed at the /i = 
value. We are here interested in values of /x/ ~ '^vr/S ~ 70 MeV, and it is a well- 
known fact that the physics at T = is independent on the chemical potential until 
a threshold value. We assume that, at yUg = and /i/ 7^ 0, this critical value is such 
that a non zero p forms with a x fixed to its /i = value. In analogy with chiral 
models, we also assume that the transition is continuous; therefore, the critical value 
/if must verify 



After some algebra, we have 



Qp2 lp=0 



(pF) = 



(4.30) 



1 

AG 



00 ^p. 



dp p 



+ 



1 



\E 







(4.31) 



vr- Jo ^ LI-^ + pFI 

Now, we consider without loss of generality positive values for /i/; we notice that 
the term (E—pi) is positive, since E is in any case greater than the constituent mass, 
which is of the order of ~ 350MeV, and we consider here ni ~ 70MeV. Therefore, 
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we can safely remove the absolute values in the latter equation. 
At the end of the calculation, we find that /if must satisfy 

As a first check, we want verify that in the chiral case, when the pion mass is 
exactly zero, the critical value /ij* vanishes. By substituting m = fif = in the 
last equation, we find back the gap equation in the chiral case. This check works 
out. Now, we want to consider the massive case; in the spirit of a l/m expansion, 
we will consider values of m such that m ^ Axo in eq. (j4.27|) . and expand all of the 
previous expressions in powers of m. At the end of this procedure, only the most 
significant term in m will be kept. We skip the intermediate calculations and show 
the final expression only 



E 



(4.32) 



I 



dp 



(p2 + A2^2F^y/2 
(4.33) 

Now, we must compute the pion mass and compare it with the latter expression. 
The static pion mass is related to the second derivative of the effective potential 
with respect to the pion field, apart a normalization factor which depends on the 
relation between canonical pion field and p field. We have 

m^=p-g^\min (4-34) 

where f-,, is the pion decay constant. 

The expression for in a NJL model without form factor and with a cut-off A is 



where M is the constituent mass. One can think of removing the cut-off and 
introducing the form factor in the expression eq. ()4.35|) by substituting M with 
(m -|- AxE"^). This procedure could appear as rather heuristic, and will be dis- 
cussed later. By doing this one finds for m^, computed at the leading order in m. 
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the following expression: 



dp 



(p2 + A2;^2^4)3/2 

At the leading order in m the equivalence 



dp p^ 



(p2 + K^xlF^f'^ 
(4.36) 



/if = mj2 (4.37) 

is therefore proved. As far as the calculation of is concerned, if we consider 
the standard NJL model with hard cut-off on the three-momentum, the described 
procedure is rigorous. In the less-used version of the model with the form factor 
some difficulties arise when computing /j^; the chosen dependence of F on the three- 
momentum only breaks covariance and a closed formula for has not been found. 
This is the reason why we did not fixed in the fit procedure of the paper |HH] , and 
why there was a little numerical mismatch between /ij* and m,r/2 in that paper. 

With a self-consistent treatment of quarks in the mean fied approximation and 
of mesons in the Random Phase Approximation, the authors of |lU3j . by using a 
NJL model with U{1)a breaking and without the form factor, have shown that the 
result of eq. ()4.37p is exact at all orders. Moreover, the result is independent on the 
size of the U{1)a breaking coupling. 



4.6 Conclusions 

In this Chapter we have continued the study of pion condensation at finite quark 
and isospin density in the NJL model that we had started in the previous Chapter 
in the case of ladder-QCD for small isospin chemical potentials. The extension 
to higher isospin chemical potentials confirms the structure predicted in ref. |77j . 
where two-fiavor QCD was simulated in the context of a random matrix model. 
Some differences between the two analyses are present, at low temperatures, in the 
region of high isospin chemical potentials, at the boundary of the region where 
color superconductivity should take place. Actually in this region we find that 
pion condensation is still active, whereas in ref. [77j the authors find that the 
pion condensate vanishes. We have also shown the expected behavior of scalar and 
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pion condensates by following different paths, for growing temperatures, both in the 
plane of quark chemical potentials /id) and in that of isospin and quark chemical 
potentials {iii,iiq). The analysis that we have performed should also be confirmed, 
with only small quantitative differences, within ladder-QCD. 
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Figure 4.1: Phase diagram for chiral symmetry restoration in the plane {piui P-d) 
of chemical potentials of quarks, at T = 0. Different regions are specified by the 
non vanishing of a given condensate, whereas the others are vanishing (p) or order 
~ m/A (xu and Xd)- Dashed lines are for the continuous vanishing of p or for 
cross-over phase transitions for Xu or Xd, whereas solid lines are for discontinuous 
behaviors. The solid lines a and b refer to specific paths at fixed values of Pq, with 
Pq = 170 MeV (line a) relative to Fig. and Pq = 210 MeV (line b) relative to 
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Figure 4.2: Critical value of the isospin chemical potential, beyond which a pseu- 
doscalar condensate forms, vs. quark chemical potential, at zero temperature. 
Hq = 204 MeV is the highest allowed value for pion condensation to occur. The 
path followed in the phase diagram of Fig. is along the upper-half of the dashed 
line in the lower half-plane. 




Figure 4.3: Scalar and pseudoscalar condensates vs. ni, for jjg = 170 MeV and 
T = 0. The path followed in the phase diagram of Fig. \4.1\ is that of the solid line 
a. 
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Figure 4.4: Scalar condensates vs. /i/ for jiq = 210 MeV and T = 0. The pseu- 
doscalar condensate is zero. The path followed in the phase diagram of Fig. \4.1\ is 
that of the solid line b. 
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Figure 4.5: Scalar and pseudoscalar condensates vs. Hd, for = 200 MeV , T = 
(see Fig. [g. 
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Figure 4.6: Phase diagram for chiral symmetry in the [fig, T) plane for zero or small 
isospin chemical potential fij. For /i/ = (central line), the cross-over transition 
line starts from the point Pt = (0, 174) and ends at the point E = (220, 85). The line 
between E and the point P^ = (293, 0) is the line for the first order transition with 
discontinuities in the (uu) and (dd) condensates. For /i/ = 30 MeV (side lines), 
the two cross-over transition lines start from the point Pt = (0, 174) and end at 
the points E^ = (190,85) and Ed = (250,85). The lines between E^ and the point 
= (263, 0) and between E^ and the point Pj^ = (323, 0) are the lines for the first 
order transitions with discontinuities in the (uu) and (dd) condensates respectively. 
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Figure 4.7: Phase diagram for chiral symmetry restoration in the plane {pu, f^d) 
of quark chemical potentials, at T = 60 MeV , which is below the temperature of 
the critical ending point (see Fig. \4-(^ - Different regions are specified by the non 
vanishing of a given condensate, whereas the others are vanishing (p) or m/ A (xu 
and Xd)- Dashed lines are lines for the continuous vanishing of p or for cross-over 
phase transitions for Xu or Xd, whereas solid lines are for discontinuous behaviors. 
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Figure 4.8: Phase diagram for chiral symmetry restoration in the plane [fiu, f^d) of 
chemical potentials of quarks, at T = 140 MeV , which is above the temperature of 
the critical ending point (see Fig. |^.6| j. Different regions are specified by the non 
vanishing of a given condensate, whereas the others are vanishing (p) or m/A (xu 
and Xd)- Dashed lines are lines for the continuous vanishing of p or for cross-over 
phase transitions for Xu or Xd, whereas solid lines are for discontinuous behaviors. 
The solid line a refer to the path at Hq = vs. fii > followed in Fig. \4-lC\ 
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Figure 4.9: Scalar and pseudoscalar condensates vs. Hd, for yU„ = 200 MeV, T = 
60 MeV (see Fig. Ul\ ). 




Figure 4.10: Scalar and pseudoscalar condensates vs. /i/, for /ig = and T = 
140 MeV . The figure is obtained following the path a in the phase diagram of Fig. 
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Figure 4.11: Scalar and pseudoscalar condensates vs. Hd, for /i„ = 200 MeV, T = 
140 MeV (see Fig. CT . 




Figure 4.12: Critical value of the isospin chemical potential, beyond which a pseu- 
doscalar condensate forms, vs. temperature, at fig = 0. 
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Figure 4.13: Pion condensates vs. T , for fig = and fij = 200 MeV . The scalar 
condensates Xu CLnd Xd O'l"^ of order ~ m/A. 
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Chapter 5 



Pion and kaon condensation in a 
3-flavor NJL model 

We analyze the phase diagram of a three-flavor Nambu-Jona-Lasinio model at finite 
temperature T and chemical potentials fid, fJ^s- We study the competition of pion 
and kaon condensation and we propose a physical situation where kaon condensation 
could be led by light quark finite densities only. We refer to the paper |lU4j . 

5.1 Introduction 

In the Chapters III and IV we have studied pion condensation by using two effective 
models, in order to mimic 2- flavors QCD. The existence of phases with condensed 
pion is confirmed by lattice analyses [7^ llU5j . within chiral models |H IH^ ITUH| llU7j 
and within microscopic models [771 EHl EZl 1103] . 

In this Chapter we extend our previous analyses to the strange flavor. It has 
been suggested a kaon condensed phase in the context of high density nuclear matter 
[SSmnHllinnillini and, more recently, in particular in the CFL regime PTIITTTIITT^ : 
in this context, the mass of the strange quark plays the role of distinguishing the 
pattern of symmetry breaking between a two flavor superconducting (2SC) and a 
three flavor color flavor locking (CFL) phase. For this reason, an exhaustive study of 
these phenomena must take into account the strange quark sector. It is also worth to 
mention that taking into account weak equilibrium and color and electric neutrality 
many other interesting superconducting phases could exist [201 [73 1113| I114j . 

These study could regard the physics of compact stars, more than that of col- 
liders. The existence of stars whose core is so dense that quarks can be treated 
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as fundamental degrees of freedom (i.e. quasi-particles) has been considered for in- 
stance in [7U I115t I116t 1117] . Actually, compact stars can be thought as the only 
hypothetical experimental settings for color superconductivity to occur. 

We will propose a physical scenario in which kaon condensation could be driven 
solely by a light quark finite density. The possibility of a kaon condensation at 
/is = could be important in the context of neutron stars, where, by neglecting 
electroweak effects, the densities are associated with light quarks only. 

This has been the first study in literature dealing with kaon condensation in a 
microscopic model^. 



5.2 The model 

Our purpose is to explore the general structure of the phase diagram for chiral 
symmetry and pion and/or kaon condensation in three-flavor QCD at non zero 
quark densities and its evolution in temperature, by using the Nambu-Jona Lasinio 
model (NJL). In this way we generalize, to three flavors, previous works concerning 
the simpler two- flavor case [7T| |HH|- To our knowledge, so far the only study of 
meson condensation in the general case of three-flavor QCD is based on a chiral 
Lagrangian |1()8^ IHIlj : nevertheless the use of a model such as the NJL model 
or ladder-QCD jH2] is necessary if we want to include the effects of a net baryon 
charge and also study chiral symmetry restoration through the behaviour of scalar 
condensates, besides the pseudoscalar ones. 

Let us thus start with the Lagrangian of the NJL model with three flavors u, d, s, 
with current masses rriu = ma = m and and chemical potentials fiu, fJ'd, fJ's 
respectively 



■£ — -^0 + + ^fi + ^int 



md'if - + + GJ2 [(^^a^)^ + (^i75A„^)^j (5.1) 



a=0 



^In Ref. jll8) a study of meson properties at finite quark densities within the NJL model is 
performed. However, they do not deal with meson condensation. 
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where 



/ u\ / fin \ / m \ 

^ = I d \ , A = I fid \ , M = I m \ 
\s J \ fis J \00m, / 



M is the quark current mass matrix which is taken diagonal and A is the matrix of 

[2 

the quark chemical potentials. As usual Xq = J - 1 and , a = 1, ...,8 are the 
Gell-Mann matrices. 

For M = and A = ot A oc I, the Lagrangian is U{3)l ® U{3)r invariant. The 
chiral symmetry is broken by M 7^ which also breaks SU{3)v down to SU{2)v as 
m 7^ rris. However, this symmetry is also lost, since we generally consider 7^ fid- 
The remaining symmetry is thus the product of three independent phase transfor- 
mations Uu{l) ® Ud{l) ® Us{l)- For simplicity, in this application we do not consider 
the 't Hooft determinant, which explicitly breaks U{1)a- Also, we turn off the 
electroweak effects and we do not consider the possibility of di-quark condensation; 
therefore, our results will be reliable only in the regime of intermediate densities 
(with chemical potentials lower than 400 ^ 500 MeV). 

We note that we can express £^ either by using the variables fiu, fJ'd, fJ^s or, by 
introducing an averaged light quark chemical potential 

fig = {fiu+ fJ'd)/'^ (5.2) 

and the three combinations 



/i= (/i„ + yUd + /is)/3; fii = {flu- fid)/'^; fiy = ifJ'q - fJ's)/2 (5.3) 

which couple to charge densities proportional to the baryon number, to the third 
component of isospin and to hypercharge respectively, 

2 

C^ = fx + fii ^Us* + ^ /iy ^Ug^ (5.4) 

Otherwise we could decide to express in terms of fig, fii, fis, where 3/ig would 
couple to the baryon number and fis would be the chemical potential associated 
with strangeness. Obviously, the two choices are equally legitimate, and we will 
adopt 
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To study whether a pion and/or kaon condensate is formed, we need to calculate 
the effective potential. To do this we generalize to SU (3) p the bosonization proce- 
dure exposed in the previous Chapter. The one-loop effective potential we get is: 



V = ^ixl + Xl + Xl + 2pl, + 2pl + 2pl) + ^iQg (5.5) 



1 """^ r d?v I -IbF^iP) A pud -i^F^{p) A 

^log = -^y, TT^ log ^ K -l^F\p) A p„ 

P W \ ^^F\p) A p., ^,F\p) A p,, 

(5.6) 



P/)7o - p-7 - ("^/ + ^^(P) A Xf) 



A^ 

where the form factor -F(p^) = — (A is a mass scale) has been coupled to 

A^ _l_ pi 

fermionic legs to mimic asymptotic freedom In eq. ()5.(j|l tr means sum over 
Dirac, flavor and color indices and ujn are the Matsubara frequencies. The dimen- 
sionless fields Xf ^-nd p^^i are connected to the condensates by the following relations 



X7 = - 4G <*'*^> 



A 

(5.7) 



Pff = - 2G 



and are variationally determined at the absolute minimum of the effective poten- 
tial. We expect a superfluid phase with condensed pions when the isospin chemical 
potential exceeds a critical value pf (pf = at T = pg = 0); analogously a 

kaon condensation phase is expected when px is higher than the kaon mass mx- 
Since we are neglecting electroweak effects, the four kaons have all the same mass 
rriK- We can express the kaon chemical potentials as a function of quark chemical 
potentials, by studying the effective action in the kaon fields sector. In practice, we 
have to study the effective kaon masses when chemical potentials are included. In 
the normal phase (i.e. without meson condensates) we find: 

Pk+ = 2py + Pi ; Pro = 2pY - Pi (5.8) 
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For /i/ = 0, at T = the critical value for kaon condensation should be /iy = mx/2. 

To fix the free parameters of the model, which are A, the average light quark 
and the strange quark masses m, rris and the coupling G = g/A"^, we calculate the 
charged pion and kaon masses (by fixing their decay constants to the experimental 
values) and the light quark scalar condensate, in the vacuum. The fit procedure is 
analogous to that of Chapter IV. 

By choosing the free parameters as follows 

A = 1000 MeV; g = G A'^ = 6; m = 1.7 MeV; m, = 42 MeV (5.9) 

we get for the light quark condensate, the constituent light quark mass, the pion mass 
and the kaon mass respectively the results (^/^/) = -(252 MeV)^ Mq = 385 MeV, 
= 142 MeV, = 500 MeV. Furthermore, we find the expected agreement 
between m,r/2 and fif (within 1%), and also between /ly and rriK/'^ (within 4%). 

The very low value we find from our fit for the quark masses is entirely due to 
the choice of form factors (instantaneous approximation), see |119j . but their ratio 
ms/m ^ 25 is in agreement with the value given by the ratio of the pseudoscalar 
masses. 

5.3 Phase diagram for chiral symmetry and me- 
son condensation 

As usual, we determine the phases of the model by looking at the absolute minimum 
of the effective potential, given in Eqs. (|5.5|) . (j5.6p . 

The effective potential depends on three scalar fields and three pseudo-scalar 
fields; moreover, a non zero value of a field pff induces an effective coupling be- 
tween the flavors //'. Therefore, in general, the minimum has to be searched in a 
six-dimensional space, and the analysis turns out to be very hard from the numerical 
side. 

Since there are four thermodynamic parameters, namely T, /i^, /is (or T, fi, /i/, /iy), 
we start from T = and then follow the evolution in temperature of two-dimensional 
slices in the space of chemical potentials. We recall that the global symmetry in the 
model is given by the product of three U{1) groups. The possible formation of a 
non-zero value of one of the pseudoscalar flelds p^^' implies the spontaneous break- 
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ing of a U{1) symmetry with the appearance of a Goldstone boson. The transition 
to this superfluid phase leaves two U{1) groups unbroken and between them there 
is always the U{l)y. For instance, the formation of a pud condensate implies the 
breaking of the U{1)y symmetry, whereas U{1)y ® U{1)y are left unbroken. At the 
threshold one of the charged pions, depending on the sign of /x/, is the Goldstone 
boson, whereas the other mesons of the octet remain massive (see also pUj). 

Thus we characterize different regions of pion or kaon condensation with the 
pseudoscalar field which acquires a non vanishing vacuum expectation value (v.e.v.), 
whereas the other v.e.v. 's are vanishing. Then, to facilitate the physical interpreta- 
tion and to directly compare it with the results of ref. jHOj, in the pictures we label 
these regions directly with the symbol of the particle that condenses, rather than 
with that of the corresponding field Pff- 

The scalar fields X/ do not break any symmetry as chiral symmetry is explicitly 
broken by the current quark masses. However these fields undergo cross-over or 
discontinuous transitions for given values of the thermodynamic parameters. 
Thus we distinguish between regions where a non zero v.e.v. of a scalar field is pri- 
marily due to dynamical effects, and we indicate them with the relative symbol x/? 
from regions where their values are of order ~ rrif/A, where we put no symbol. We 
also avoid putting the symbol x/ in a region where the pseudoscalar field associated 
with that flavor starts to form, and Xf is not yet of order ~ rrif/A. Furthermore, 
both for pseudoscalar meson condensation and for the transitions associated with 
scalar fields, the boundaries of the regions separated by discontinuous transitions 
are represented by solid lines whereas dashed lines indicate continuous transitions. 
To start with, we set /!„ = —pd, which is the case studied by Toublan and Kogut 
in ref. [00] • However, to compare our results with theirs, we have to remark that 
we adopt different combinations of the microscopic chemical potentials /i^. 
In this case, where Pq = (see Eq. ()5.2p ). we plot the phase diagram in the plane 
of fiy and fij, which are proportional respectively to ps and to pu — Pd and which 
correspond to half of the chemical potentials employed in ref. pnj . 
In Fig. 15.11 we plot the quadrant with /iy > 0,yU/ > of the phase diagram at zero 
temperature. Let us start by looking at what happens by moving along the vertical 
line at /i/ = 0^, where O"*" stands for an infinitesimally small positive value, in order 
to remove the degeneracy between and K^. 
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In this case one starts from a "normal" phase (see also ref. [HO]) characterized by 
chiral symmetry breaking by large scalar condensates and then, for yUy ~ ~ 
250 MeV, the scalar Xu and Xs condensates smoothly rotate into a p^g condensate, 
whereas the Xd fi^ld remains unchanged across the transition. We have indicated in 
the picture this new region with the symbol {K^) (rather than with p„s) as with 
these signs of the chemical potentials it is the positively charged kaon to condense. 
By further increasing [ly the dynamical effect associated with the strange quark 
condensation is weakened by a high strangeness density and the favored phase be- 
comes that with Xs ~ ^s/^ and large XuiXd- The kaon condensate vanishes in 
this case through a first order transition. The same situation hereby described also 
holds when we move along vertical lines with pi < m7r/2, with a linear decrease 
of the critical /ly for kaon condensation and a weak linear increase for the second 
transition at higher /xy. 

On the other hand, when we move along the horizontal line at /iy = and we 
increase /i/, we pass from the "normal" phase to a pion condensation phase (in this 
case the 7r+) described by a finite pud-, whereas Xs remains large. This situation 
was already known by the previous studies of the two- flavor NJL model [HE]- This 
transition is also second order. This situation persists at larger pi too. In the in- 
termediate regions there is a competition between pion and kaon condensation. For 
instance, by fixing pi > starting from low py in the 7r+-condensed phase with 

a large Xs and increasing /iy^, we find that the system goes to the phase where the 
condenses when py > m/f/2. Consequently Xs becomes ~ rris/A whereas Xd 
becomes large (if pi < 325 MeV) and the full transition is discontinuous. When py 
is about 300 MeV (with a linear weak increase for growing pj), the kaon conden- 
sate vanishes and a pion condensed phase shows up, with all the scalar condensates 
Xf ~ rrif/A. The transition is first order. Finally, for high pi and py in the range 
of values corresponding to kaon condensation, there is a phase analogous to the 
previous one, with a kaon condensate instead of a pion condensate and Xf ~ rrif/A. 
This region is connected to the others through discontinuous transitions too. 

Along the line pi = the Nielsen-Chadha theorem [HHI would apply. The con- 
densation of a kaon field would cause the following pattern of symmetry breaking 



^In this case, working at /iq — 0, it simply means to increase |^s|. 
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m 



SU{2)0U{1) ^U{1) (5.10) 

with three broken generators. Due to their dispersion relations, there are only 
two corresponding Goldstone bosons 

In Fig. 15.21 we plot the behaviour of the scalar and pseudoscalar condensates 
vs. fij at fiq = 0, /iy = 230 MeV and T = corresponding in Fig. 15. II to the path 
marked by the solid line a. It is worth remarking the rotation of the Xu and Xs fields 
into a Pus field at the continuous transition for kaon condensation (at fij ^ 50 MeV) 
and the further "exchange" of the two pseudoscalar fields pus and pud at the pion 
condensation transition (at pi ~ 125 MeV). 

In Fig. 15.31 we plot the behaviour of the scalar and pseudoscalar condensates vs. 
Py Sit pq = 0, Pi = 200 MeV and T = (following the path of the solid line b in 
Fig. EH)- 

In Fig. 15.41 the phase diagram of Fig. 15.11 is extended to negative values of pi 
and py too. This picture is thus simply obtained by refiection of Fig. l5.1l aroimd its 
axes at Py = and pj = 0. The only difference is that in the other three quadrants 
different pseudoscalar mesons of the octet condense (see also ref. pU]). 

The evolution of the scalar fields for growing temperatures is well-known. In 
particular, as expected the growth of temperature fights against pion and kaon 
condensation. We find that the kaon condensed phase disappears at /i/ = for 
T > 85 MeV. In Fig. 15.51 we show the evolution of the phase diagram of Fig. 15.11 
at T = 100 MeV. The region for kaon condensation has shrunk and the transitions 
associated with chiral symmetry approximate restoration have become continuous. 
Above T ~ 110 MeV the regions of kaon condensation disappear. In Fig. 15. 61 we plot 
the further evolution of the phase diagram of Fig. 15. H at T = 140 MeV. In the range 
of values of pi,pY that we have considered, only the regions of pion condensation 
and the usual regions with chiral symmetry breaking separated by continuous tran- 
sitions from analogous phases with Xs ~ fng/A for large py remain. Finally, pion 
condensation disappears above the cross-over critical temperature T = 192 MeV. 
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5.4 Kaon condensation in the region of high /ig 



So far we have considered the hmit /i„ = —fid (and thus fig = 0). A more general 
picture of the full phase diagram can be grasped by considering two other cases, 
which are fiu = fid (and thus fij = 0), and fis = 0, which has been already considered 
in ref. jHHl- In this section, we suggest a hypothesis for a kaon condensation driven 
mainly by light quark finite densities, with fis low or even zero. We recall that in the 
case of /is = there is an equal number of strange quarks and anti-quarks, and in 
that case kaon condensation would concern a light particle (antiparticle), associated 
to the external field (the chemical potential), and a strange anti-quark (quark) 
belonging to the sea. In ref. |HH]; we had not inserted the pseudoscalar fields with 
strangeness content: the possibility of kaon condensation at = would modify 
the phase diagram in the {fiu, f^d) plane of ref. |HH1 with the insertion of regions 
with condensed kaons. However, except for this, no further modification would be 
necessary, whereas the possibility of kaon condensation at zero or low values of fig 
can be also studied for fij = 0, which is therefore the last situation we examine. 
From its definition in Eq. ()5.3|) . we see that fiy measures the difference between 
the strange quark chemical potential and fig (which at /i/ = is the chemical 
potential associated with one of the light flavors) similarly to fij which measures 
the unbalance between the u and d quark chemical potentials. Therefore, since pion 
condensation occurs, at T = and for instance at fid = 0, when fiu > m^r, similarly 
we could expect a kaon condensed phase at T = and fig = when fig > rriK- 
However, within our approximations, this does not happen, the reason being that 
riix is higher than the critical value of fig for the melting of the dynamical part of 
{^f'^f), f = u,d. In other words, when the condensates of the light quarks start 
to decrease to values ~ m/A, the effective kaon mass starts to increase and kaon 
condensation always remains unfavored (the same thing does not happen at /i^ = 
because the melting of the dynamical part of (ss) occurs at higher values of fig and 
actually we find kaon condensation as shown in the last section). Nevertheless, since 
we do not have taken into account the effects of di-quark condensates which could 
break chiral symmetry (as it happens for instance in the CFL phase), it would be 
interesting to further analyze what happens with their inclusion. We will analyze 
this case in further work. 
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For the time being, let us show, in Fig. 15.71 the T = phase diagram in the 
[fig, fig) plane at fij = 0~, where 0~ stands for an infinitesimally small negative 
value, which is necessary to select which are the kaons that condense (in any case, 
within a neutron star a finite negative isospin chemical potential is present). The 
structure of this phase diagrams recalls that of [HH] for pion condensation in the 
[fJ-u, fJ-d) plane, with a "normal" phase in the middle, characterized by high values of 
the scalar condensates Xu, Xd, Xs- This phase is separated from the regions of chiral 
symmetry approximate restoration (relative to one or more fiavors depending on the 
values of chemical potentials) which are characterized by the respective scalar field 
Xf ~ rrif/A. Furthermore this region is separated through discontinuous transitions, 
both from the "normal" phase and from the regions of kaon condensation. There are 
four different regions of this kind, distinguishable by the type of the condensing kaon 
and by the values of the associated scalar condensates. Finally, the "normal" phase 
turns into a region of kaon condensation through the smooth rotation of Xs and one 
of the light flavors (according to the sign of fiy one between Xu or Xd remains large) 
into a kaon condensate, when /iy > mK/2, with a second order phase transition. 
We can see that the presence of a high value of lowers the corresponding \fis\ 
necessary to have kaon condensation: moreover, a value of ;U/ 7^ would favor a 
kaon condensation with lower values of \fiq\ and/or \fis\. 

Since we do not know exactly the critical value fi^ for chiral symmetry restoration, 
and since we are interested here in the case where fig > mx, instead of varying /i^ 
(for instance performing a different flt) we take the kaon mass as a free parameter, 
and we present, in Fig. 15.81 the same result as in Fig. 15.71 for ttlk = 300 MeV, 
which is the upper limit for having a kaon condensation at Hs = 0. 



5.5 Conclusions and further developments 

In this Chapter we have presented a calculation of the QCD phase diagram obtained 
by using a three-flavor NJL model, without U{1)a breaking. We have considered the 
plane (/i/, /iy) to study the competition of kaon and pion condensation at zero and 
flnite temperature and to establish a comparison with a previous analysis performed 
in a chiral model The results of the microscopic model accurately reproduce 
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those achieved in the chiral model only for low densities, since in the high strangeness 
density regime approximate chiral symmetry restoration (relative to the strange 
sector) occurs, disfavoring kaon with respect to pion condensation. Moreover we have 
suggested the possibility of kaon condensation driven only by light quark densities: 
this situation could be of interest for the core of compact stars, where the densities 
could be approximatively associated with light flavors only. 

To make this study more quantitative, we should include the effects we have up 
to now neglected. In the first place, we have to consider di-quark condensation in 
the regime of high densities; those studies have been driven in recent years in the 
context of microscopic models. The last goal of this project would be taking into 
account meson condensation in superconducting regimes. The possibility of kaon 
condensation in the CFL phase (KO-CFL) has been recently proposed in the 
framework of a chiral model. It would be interesting to perform the same study 
but with a more fundamental model: the problem is that the structure of the one 
loop term of the effective potential results rather hard to deal with, even though 
conceptually simple. Actually, we have to compute the determinant of a 72x72 
matrix, corresponding to 3 colors, 3 flavors, 4 Dirac, 2 Nambu-Gorkov indices. In 
the Nambu-Gorkov formalism one introduces a bi-spinor 

The result is that di-quarks can be included in the path-integral formalism by means 
of the calculation of a fermionic integral in Nambu-Gorkov space. The price to pay 
is the doubled dimension of the spinor itself. 

For di-quark fields, we should consider the following combinations: 

■522; Sss, S77 (5.12) 

with 

SAA' = (^"^CtsTaAa'^) (5.13) 

where r^, A^' are Gell-Mann matrices acting in flavor, color spaces respectively. In 
the pure CFL regime, these fields are all equal; here we consider that, due to chem- 
ical potential asymmetries or to the condensation of a flavor-mixing pseudoscalar 
field, they can be different. In order to get the ground state, we have to minimize 
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the effective potential with respect to at least seven variables (3 scalar condensates, 
3 di-quark condensates, 1 meson condensate). Actually, pseudoscalar and di-quark 
condensates mix different ffavors and colours, which turn out to be variables depen- 
dent on each other. Obviously, the study gets cumbersome from the numerical side. 
Once overcome this hurdle, it would be possible to study whether a kaon conden- 
sate could arise, at high densities, by considering only light quark densities. This 
would certainly be a remarkable results for the physics of stars. The competition 
of di-quark and meson condensates in a NJL model has been partially afforded in 
the papers jl2()[ I12H I122j which followed our work. Eventually, the inclusion of 
the /3-equilibrium as well as charge and color neutrality, should be important for 
reproducing theoretically the settings attained in the core of a compact star. The 
inclusion of the 't Hooft determinant in the model would not dramatically modify 
the results we have obtained; moreover the restoration of U{1)a in the regime of 
high densities is expected. 
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Figure 5.1: Phase diagram for chiral symmetry restoration and meson condensation 
in the plane (fij, ) (it = and T = 0. Different regions are specified by the non 
vanishing of a given condensate, whereas the others are vanishing (pud,Pus) or order 
~ mf/A (Xu,Xd,Xs)- Dashed lines are for the continous vanishing of pseudoscalar 
fields, whereas solid lines are for discontinuous behaviours. The solid lines a and b 
refer to specific paths at fixed values of fiy = 230 MeV (line a) and fij = 200 MeV 
(line b). 
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Figure 5.2: Scalar and pseudoscalar condensates vs. for T = 0, fiq = and 
/iy = 230 MeV. The path followed in the phase diagram of Fig. I5.il is that of the 
solid line a. 
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T = ; f^i = 200 MeV ; fiq = 
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Figure 5.3: Scalar and pseudoscalar condensates vs. hy, for T = 0, /i, = and 
/i/ = 200 MeV . The path followed in the phase diagram of Fig. I5.il is that of the 
solid line b. 
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Figure 5.4: Phase diagram for chiral symmetry restoration and meson condensation 
in the plane (iii^^y) at fig = and T = 0. Different regions are specified by the 
non vanishing of a given condensate, whereas the others are vanishing (pud, Pus, Pds) 
or order ~ mf/A (xu, Xd, Xs)- Dashed lines are for the continous vanishing of pseu- 
doscalar fields, whereas solid lines are for discontinuous behaviours. 
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T = 100 MeV ; /ig = 




Figure 5.5: Phase diagram for chiral symmetry restoration and meson condensation 
in the plane (fii,fiY) at fig = and T = 100 MeV, which is above the temperature 
of the tricritical point. Different regions are specified by the non vanishing of a given 
condensate, whereas the others are vanishing (pud, Pus) or order ~ m//A (xu, Xd, Xs)- 
Dashed lines are for the continous vanishing of pseudoscalar fields or for cross-over 
transitions for scalar fieds, whereas solid lines are for discontinuous behaviours. 
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T = 140 MeV ; fiq = 
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Figure 5.6: Phase diagram for chiral symmetry restoration and pion condensation in 
the plane (fj,i,fj,Y) at fig = and T = 140 MeV, which is above the highest temper- 
ature ~ 110 MeV to have kaon condensation. Different regions are specified by the 
non vanishing of a given condensate, whereas the others are vanishing (pud) or order 
~ rrif/K (Xu,Xd,Xs)- Dashed lines are for the continous vanishing of pseudoscalar 
fields or for cross-over transitions for scalar fieds. 
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Figure 5.7: Phase diagram for chiral symmetry restoration and kaon condensation 
in the plane (fig, fig) at fii = 0" and T = 0. Different regions are specified by the non 
vanishing of a given condensate, whereas the others are vanishing (pus,Pds) or order 
~ rUf/A (XuiXdiXs)- Dashed lines are for the continous vanishing of pseudoscalar 
fields, whereas solid lines are for discontinuous behaviours. 
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Figure 5.8: Phase diagram for chiral symmetry restoration and kaon condensation 
in the plane (^q,^s) at [ii = 0~ and T = 0, and for rriK = 300 MeV. Different 
regions are specified by the non vanishing of a given condensate, whereas the others 
are vanishing (pus,Pds) or order ~ rrif/A (Xu,Xd,Xs)- Dashed lines are for the 
continous vanishing of pseudoscalar fields, whereas solid lines are for discontinuous 
behaviours. 
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Chapter 6 



A study of the critical line in a 
3-flavor NJL model with U{1)a 
breaking 

We employ a 3 flavor NJL model to stress some general remarks about the QCD 
critical line. The dependence of the critical curve on fig = + /id)/2 and 

= ifJ'u — /Ud)/2 is discussed. The quark masses are varied to conflrm that, in 
agreement with universality arguments, the order of transition depends on the num- 
ber of active flavors Nj. The slope of the critical curve vs. chemical potential is 
studied as a function of Nj. We compare our results with those recently obtained 
in lattice simulations to establish a comparison among different models. We refer 
to the paper |123j . 

6.1 Introduction 

In recent years, the study of the QCD phase diagram by means of numerical lattice 
simulations has improved considerably. In particular, the presence of a critical end- 
ing point, flrst discovered within microscopical effective models [21112111101, appears 
to be a solid feature pH lll24j . althought its exact location along the critical curve is 
still controversial. Moreover, the incoming realization of LHC enhances the interest 
of the scientiflc community toward this topic. 

The main problem regarding the lattice study of the QCD phase diagram is related 
to the so called sign problem; the fermionic determinant is not positive deflnite 
at flnite baryonic chemical potential, and therefore, to overcome this unwelcome 
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feature, some suitable tricks are needed. Without going into detail, the methods 
commonly used are the study of QCD at imaginary chemical potential jHSl UHl 11^5] , 
the reweighting procedure [^3] and the Taylor expansion in /i/T |^ 1126] . 

The aim of this paper is to offer an overview of results concerning the behaviour 
of the critical line, obtained in the NJL model and directly comparable with recent 
lattice analyses. In its strong simplicity, the NJL model recovers the basic structure 
of the non perturbative dynamics ruling the problem. Therefore, it can be trusted 
as a good toy model for the study of QCD phase diagram. 



6.2 The model 

Let us consider the Lagrangian of the NJL model with three flavors u,d,s, with 
current masses m„ = = m and and chemical potentials /i^, fig respectively 











d 






{ s ) 





£ — Co + Cm + Cfj^ + C4 + Cq 

8 

= md"^ + + gJ2 [(^Aa*)^ + {m-f^Xa'^Y (6-1) 

a=0 

+ [det^(l + 75)^ + detf (1 -75)^] 

where 

flu \ / m 

0/idO , M = I m \ (6.2) 
fis J \ m, 

Ai is the current quark mass matrix which is taken diagonal and A is the matrix of 

the quark chemical potentials. As usual Aq = y - I and , a = 1, ...,8 are the 
Gell-Mann matrices. 

With respect to our previous studies, we consider here the 't Hooft term which 
breaks U{1)a- It is the determinant term with coupling K in eq. 1)6.11) : for the three 
flavors case, it corresponds to a six fermion interaction. In the chiral case this term 
breaks the axial current conservation 

d^J^ = -4Nf K Im(det(^(l - 75)^)) (6.3) 
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so that, by comparing with eq. (|1.4p . we find for the topological charge density in 
this model [HI 



Q(x) = -2ir Im det(^(l -75)^)) (6.4) 

By working at the mean field level, the six fermion term can be recast into an 
effective four-fermion one. In such a way the Lagrangian reduces to the usual 
NJL Lagrangian, apart from a redefinition of the four-fermion coupling constant G 
into a new set of effective ones, taking into account the fiavor mixing arising from the 
't Hooft term 110] • Each channel in the interaction (corresponding to different 
Xa, and scalar /pseudoscalar sectors) has a different coupling, which depends on G, 
K and on the average value of other fields. 

In the following, we will work at the mean-field approximation. Therefore, be- 
cause we are now dealing with four fermion interactions only, we can calculate the 
effective potential by following a similar procedure to that of Chapter four, taking 
care of the new set of couplings. The explicit calculation will not be shown. 
If we limit ourselves to consider the three scalar condensates, and the pseudoscalar 
condensate in the light quark sector only, the one-loop effective potential we get is: 

y = ^(^l + xl + Xs) - ^^Y^ Xu Xd Xs + ^(1 - ^^^)(pL) + '^log (6-5) 



where 



log 



n=+oo 



(2^ 



K 

tr log I 75 A Pud (1 - 




KAxs 
8G2 



-75 A P«d (1 
hd 




8G2 








(6.6) 



hf = {iujn + fJ^f)^o - p-7 - Mf 

In eq. (16. 6j) tr means trace over Dirac, flavor and color indices and Un = {2n + l)7r/ (3 
are the Matsubara frequencies. The dimensionless flelds Xf and pud are connected 
to the scalar and pseudoscalar condensates respectively by the following relations 
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Pud — — 2G 



A 

(6.7) 

A 



and are variationally determined at the absolute minimum of the effective potential. 
The constituent quark masses are 

M, = m, + Ax^-A'^^ {^^3^k) (6.8) 

As an effect of the 't Hooft term, from eq. ()6.8|l we can see that the constituent 
masses receive a flavor-mixing contribution from all of the three flavors. 

For the present application, we will consider the "standard" NJL model, without 
the form factor we had considered in the two previous Chapters. Since the model 
turns out to be non-renormalizable ([G] = energy^^, [K] = energy"^), we introduce 
a hard cut-off A on the three-momentum. 

Through this Chapter, when we consider the physical case (with realistic values 
of meson masses and decay consant), we assume for the parameters the same values 
as in ref. 

A = 631.4 MeV; G A^ = 3.67; K A^ = -9.29; (6.9) 



m = — - = 5.5 MeV; m, = 135.7 MeV (6.10) 

To investigate regimes different than the real one, with a varying number of 
Nf massless flavors, we will treat the quark masses as free parameters, by keeping 
coupling and cut-off scale fixed as in eq. (j6.9|) . 

In the following, it will turn out to be more convenient to introduce the following 
linear combinations of chemical potentials: 

Aig = (/i« + /^d)/2; /i/ = (/i„ - /id)/2; (6.11) 

The quark chemical potential /ig is just one third of the baryon chemical potential 
Pq = Pb/^- 
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6.3 Behaviour of critical lines 



6.3.1 Critical temperature dependence on baryon/isospin 
chemical potentials 

The aim of this section is to shed some hght on the physics of QCD at finite baryon 
chemical potential {fig = /is/S), by comparing the physics at fig ^ and fij = 
with that at fig = and /i/ 7^ 0. In fact, the latter case can be studied on the lattice 
by means of standard importance sampling techniques. The connection of these 
two regimes could give a deeper understanding of the sign problem in the fermion 
determinant, and provide us with some procedure to check present simulations and 
possibly improve numerical algorithms. 

We will consider here the dependence of the critical temperatures 

TcifXg) = rc(/ig; ni = 0) ; Tc(/i/) = Tc(/i/; /ig = 0) (6.12) 

for low chemical potentials, obtained by a mean-field analysis of the N JL model. Ob- 
viously, mesons and baryons (and di-quarks) carry different spin and charges, and 
their properties depend differently on fig and fij] for these reasons, one could expect 
the two curves Tc{fig) and Tc{fii), when starting by the same value at fig = fij = 0, 
to be different, at least in the regime where bound states heavily influence the ther- 
modynamics of the system. On the other hand, when the free-energy is mainly ruled 
by the constituent quarks, there is no reason to expect a dependence of the critical 
curve on the sign of chemical potentials (which fixes the sign of the total charges 
associated with the system). 

Here, at the mean-field level, the effect of bound states is considered when we admit 
the formation of a pion condensate. Actually, in agreement with chiral models anal- 
yses, in the NJL model the pion effective mass dependence on chemical potentials 
can be analitically computed |lUlj (see Chapter IV). The charged pions chemical 
potential is exactly the double of the isospin chemical potential. For this reason, as 
fii is higher than some critical value at T = 0), a pion condensate starts to 

form; a similar effect happens when fig is higher than the critical value for di-quark 
condensation, which is expected to occur at values fig > 400 -r- 500 MeV (of course, 
before di-quark condensation, for /^^ ~ m^v/S and low temperatures there should 
be the liquid-gas transition for the nucleons). Since in the following we will be 
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interested mainly in the regime of relatively small chemical potentials (lower than 
~ 200 MeV) we will neglect the latter possibility. 

In fact, when the pion condensate is zero, the mean field effective potential is sym- 
metric under fiu — > — yUn, yUd — > this implies that fig ^ ±fii is a symmetry of 
the problem. Therefore, for zero p, the two curves Tc{pq) and Tc{fj,i) have the same 
analytical dependence. 

In Fig. 16.11 16.21 the phase diagrams in (/iq,T), {fii,T) spaces are shown; starting 
from the common value Tq = 201 MeV, corresponding to fig = fij = 0, the cross-over 
curves coincide up to the value of about 150 MeV for both chemical potentials. For 
higher /i/ (and temperatures lower than ~ 200 MeV) we are in the condensed pions 
phase; in agreement with |101j this regime will persist until fij ~ 860 MeV, before 
the saturation regime takes place. In that regime the pion condensate vanishes. The 
line Tc(yU/) which defines the pion condensation phase should be in agreement with 
the calculation coming from a relativistic, interacting Bose gas. This line should be 
continuously connected with the critical curve of a BCS system, which should take 
place at higher chemical potential, according to perturbative QCD. The cross-over 
transition, according to the papers of Son and Stephanov [HEII, should take place 
at a value /i/ ~ 400MeV. The saturation regime (i.e. the regime at high fij where 
the pion condensate vanishes), which is not expected from a perturbative calcula- 
tion in the asymptotic high /i/ regime, is a feature of the present model, which is 
not reliable at high However, it is a feature common to the Random Matrix 
approach [77j . 

On the other hand, if we follow the cross-over line Tc{fig), we will find, as ex- 
pected, a critical ending point for /ig = 330 MeV, T = 42 MeV, and a line of first 
order transitions for higher fig. 

By summarizing, at this level of calculation (mean field) , as long as isospin chemical 
potential is lower than the critical value for pion condensation the curves Tc{fiq) and 
Tc{fii) have the same analytical expression Tc{fig) = Tc{fii). This could be interest- 
ing in the attempt to extend lattice results from pi ^ and /Xg = to /Xg 7^ and 
/i/ = (at least in the region of low chemical potentials). Our conclusions appear 
to be in agreement with the authors of ref. |128j (both from the lattice and from a 
hadron resonances gas model) . 
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Figure 6.1: Phase diagram in the plane (fig,T), for the physical 2+1 case; /i/ and 
fig di^e set to zero. Dashed/solid lines indicate cross-over/first order transitions; 
consequently, the dot in the picture labels the critical ending point. 
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Figure 6.2: Phase diagram in the plane {fif,T), for the physical 2+1 case; /ig and 
fig are set to zero. The dashed line indicates a cross-over transition for the scalar 
condensates, whereas dotted line stands for genuine second order transition for the 
pion condensate. In the regions of a non-vanishing pion condensate, the discontin- 
uous behaviour of the scalar condensate turns into a continuous one, therefore the 
critical ending point is not present in this case. 



6.3.2 Order of the transition by varying rus 

It is generally assumed from universality arguments [2Zj and lattice analyses PH] . 
that the order of the phase transition by increasing the temperature at zero density 
can change as the quark mass values are varied. If in the realistic case (with physical 
quark masses) the zero density transition is expected to be a cross-over, lattice 
analyses seem to show a first order transition when the three light quark masses 
are small enough. In particular, by taking = = 0, there should be a critical 
value for below which the transition turns into a discontinuous one: different 
lattice approaches find to be half of the physical value of the strange quark 
mass |129t ll3Uj or ~ 5 ^ 10 niu^ (physical) 
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To study this aspect in the NJL model, we start from the parameters fit of [T^ . 
and we take the quark masses as free parameters; namely, we take the four and 
six fermion couplings fixed so as to reproduce the phenomenology of the realistical 
physical situation. This way of proceeding could seem as rather arbitrary, but for any 
value of the masses considered here, we have verified that the output parameters 
have reasonable values (critical temperatures 130 -r 200 MeV, light quark scalar 
condensates (—250 -. — 240 MeV)^, constituent light quark masses 250 -i- 350 MeV). 
Of course our results will be strongly model dependent, both for the choice of a 
specific set of parameters and of a particular model itself. For instance, the NJL 
model provides an estimate of the position of the critical point at lower temperatures 
and higher chemical potentials with respect to those obtained in ladder-QCD and 
the ones from recent lattice simulations |124j . 

We can clearly understand the appearance of a first-order phase transition along 
the temperature axis if we consider a Ginzburg-Landau expansion of the free energy 

F($) = + 6$^ + - m$ (6.13) 

we are here considering three fiavors with the same mass m, and $ stands for 
the common value of the scalar fields. The relevant feature of the instanton-induced 
coupling of eq. ()6.5|) is the appearance of a cubic term in the expansion. At T = /i = 
we have a < 0, c < 0. Following universality arguments jSI], this can lead the 
transition to be first-order in the chiral case. This happens if the cubic term is still 
negative at the temperature where a changes sign and becomes positive. A finite 
mass term plays the role of "smoothing" the cubic term, so that for m > the 
transition turns into a cross-over. This study will be the subject of a future work. 

In Fig. 16.31 we show the behaviour of the light quark scalar condensates vs. 
temperature, in the m„ = = limit, and for vanishing chemical potentials too. 
In the upper picture, is taken to be zero, and for a temperature of about 130 MeV 
there is a sharp first order transition. As we increase rris the discontinuity of the 
scalar condensates reduces (in the lower picture the case = 8 MeV is shown), and 
when rris exceeds the critical value = 10 MeV, the zero temperature transition 
turns into a genuine second order transition. The value we get for is in any case 
smaller than the one from lattice predictions. 

In Fig. 16.41 we plot the phase diagram in the {fig, rris) space, by taking m„ = 
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fixed to zero, and fij = fis = 0. The label //// indicates, for every couple {^q,ms), 
whether, by increasing the temperature starting from zero, the transition is a first 
or a second order one. Actually, up to rris < m^' = 10 MeV, we have first order 
transitions for every value of fig, and consequently there is no critical point; for rris 
slightly above the critical value, the critical point locates at a value of fig ~ 200 MeV. 
The critical value for fig grows together with rris until the strange quark decouples 
from the two light quarks and the critical point /ig coordinate is independent on rris, 
and lies at fig ~ 300 MeV. 

Finally we have studied whether, in agreement with lattice analyses, a first order 
transition persists when we consider a non zero but small m„ = md- this does 
not happen in the NJL model with our choice of parameters, for any value of rris. 
A recent work based on the linear sigma model has found the critical value for 
rriu = rrid = rris to be merit = 40 ± 20 MeV |131j . 



6.3.3 Critical lines as a function of Nf 

Even though lattice analyses at finite density still present ambiguities in their dif- 
ferent approaches, some general features about the critical line appear to be rather 
solid. In particular, if Tq is the critical temperature for zero chemical potentials, 
the dependence of (T/Tq) as a function of (/u/Tq) should be parabolic, (at least in 
the regime (/i/To) < 1), i.e. of the form (T/Tq) = 1 — a^fi/Toy. Secondly, the a 
coefficient should depend on the number of flavors Nf, increasing with Nf] in fact, 
the curves relative to Nf = 2, 2 + 1, 3 should be very close to each other and the one 
relative to Nf = 4 should be steeper. 

It is clear that a dependence of the a coefficient on Nf must be related in any case 
to a coupling between the flavors; otherwise, the effective potential would become a 
sum of single flavors contributions, and the critical temperature would not depend 
on Nf. This fact can give us an idea about the strength of the coupling between 
the flavors at the phase transition, and of its possible reduction with temperature. 
We will check the issue relative to the cases 2, 2 + 1, 3 in the framework of the NJL 
model. 
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Figure 6.3: Behaviour of the light quark scalar condensates as a function of temper- 
ature at zero chemical potentials and m„ = md = 0, for = (upper picture) and 
nig = 8 MeV (lower picture). In the upper picture, x = Xu = Xd = Xs, in the lower 
X = Xu = Xd- When the strange quark mass exceeds the critical value = 10 MeV, 
the discontinuous behaviour turns into a genuine phase transition. 
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Figure 6.4: Phase diagram in the (fiq,ms) space, at = nid = and fii = fig = 
0. For every region in the diagram, the label I /II means whether, by increasing 
temperature and starting from T = 0, the transition is of first or second order. 
Actually, the line in the diagram separating the two different regions follows the 
critical point by varying . 
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In this context, the aforementioned situations labeled by 2,3 are those with 2,3 mass- 
less flavors (in the case 2, is set to 5 GeV to decouple the strange quark); 2+1 
is the physical case with realistic values of quark masses. In the following, /i will 
indicate a common value for the chemical potential equal for all the active flavors. 



T/To 

1 
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0.5 1 

Figure 6.5: Plot of the three cases 2,2 + 1,3 in dimensionless units (ii/Tq,T /Tq), 
for yu/To < 1. The curves relative to 2 and 2 + 1 (the dashed one) almost overlap. 
In disagreement with lattice simulations, the curve relative to 3 stands slight above 
the others. 

In Fig. 16. 51 we show the phase diagrams relative to the cases 2 + 1, 2, 3. Obviously, 
in the case 3, we are in the situation where < and we have only first order 
transitions. As we are varying the quark masses to consider different situations, Tq 
has a large range of variation (from 130 to 200 MeV); it is impressive that, as we 
plot the phase diagrams in dimensionless units, these large differences almost cancel 
out. This is the most striking evidence that our approach, in its simplicity, has some 
validity. 

When showing the three cases on the same diagram, we can clearly observe the 
surprising overlap of the results. The agreement between 2 and 2 + 1 cases (apart 
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from numerical instabilities related with 2 + 1 case) is remarkable, and good between 
2 and 3 cases; the slope of case 3, with respect to the lattice previsions, is smaller 
than that of 2, but the difference is slight. 



The results we get for a are the following: 



a = 0.1995 ±0.025 



2 + 1 a = 0.2496 ± 0.0623 (6.14) 



3 a = 0.1614 + 0.011 

The worse precision we obtain for the 2 + 1 case depends on the larger error in 
the determination of the cross-over curve. 

For completeness, we quote the results for a from lattice analyses: 
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We find that our predictions for a are bigger than those obtained by lattice 
approaches, apart from |126j for A^j = 3; in this case, the results are comparable. A 
recent study based on a hadron resonance gas model |128| gives the result 
a = 0.17 + 0.01 and this value is in a good agreement with our results. A study 
within a chiral quark model gives for a a value of about 0.1, extracted from Fig. 2 
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of [HI]. A Statistical Hadronization model provides the result a = 0.2299 jl32j : in 
this case, the critical line must be associated with the freeze-out of hadrons rather 
than with chiral symmetry restoration. 

We have also studied the behaviour of the critical curve in ladder-QCD [82 (in 
its version [7T]); in this model there is no coupling between flavors, therefore it is 
independent on Nf. We have found that the critical curve is flatter than that of 
NJL model, namely the coefficient a is much smaller, and hence closer to lattice 
predictions: a = 0.0797 ± 0.0056. We have attributed this feature to the lack 
of coupling between flavors in the model. On the other hand, since the value of 
the a coefficient that we find in the NJL model is slightly higher than the value 
found in ref. |128j and sensibly higher than the values obtained from other lattice 
analyses, we can argue that by introducing an effective reduction of the K coupling 
with temperature, as considered in ref. [T^ll33j . a would decrease at the same time. 
Therefore, by considering the following temperature dependence of the 't Hooft term 

K{T) = Ko exp(-T/Ti)2 (6.16) 

we have verified that reducing Ti the critical temperature at /i = is also reduced; in 
this way, the curve gets fiatter. By taking two different values for Ti, and considering 
for simplicity the case 2, we find the value a = 0.1186 ± 0.0061 in the case 
Ti = 160 MeV, and a = 0.1084 ± 0.0035 in the case Ti = 100 MeV. In this way, the 
agreement with lattice is considerably improved; this can be considered an indirect 
proof of U{1)a effective restoration with temperature. 

For the sake of completeness, we have also studied the two flavor NJL model with 
the 't Hooft determinant. The model is not completely equivalent to the three flavor 
model in the limit of inflnite rus, since in the latter case a strange quark loop gives 
a contribution, proportional to the strange quark condensate, to the light quark 
constituent mass. In any case we do not expect a dramatic change in our results, 
with respect to previous expectations, for the behaviour of the critical curve. 
We consider in this case the following expression for the interaction part of the 
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Lagrangian: 



Ant = G'i[(gg)^+(gfg)2+(gi75g)2+(gi75fg)2]+G2[(gg)^-(gTg)^-(gi75g)2+(gi75rg)2] 

(6.17) 

with 

Gi = {l-p)Go ■ G2 = PGo (6.18) 

For the two flavor case, the 't Hooft term writes as a four fermion interaction 
with coupling G2. By looking at the G2 term of eq. ()6.17|l . we see that combinations 
like 

{qqf - {qiim? (6.19) 

are not invariant under U{1)a transformations. On the other hand the Gi term is 
SU{2)l SU{2)fi invariant. The /3 coefficient tells us how hard the flavor mixing 
is; it is maximal for Gi = 0, namely for /3 = 1. 

For the choice of the parameters Gi and G2 we follow the approach proposed 
in ref . fH^ |HZj . In ref . ^] the authors study the original two flavor Lagrangian, 
proposed by Nambu and Jona-Lasinio with Gi = G2 and therefore /3 = 0.5. 
The value we find in this case for a is very similar with the result we obtained in 
the SU{3) case in the limit —>■ 00: 

a = 0.2107 ±0.0214 (6.20) 

The authors of ref. jHZI take instead /5 as a free parameter. Here we furthermore 
consider the possible dependence of the G2 coefficient on the temperature, 
G2 = G2{T = 0) exp(-(T/Ti)2). 

If we take G2 independent on the temperature (namely with Ti = 00), the value 
we find for the a coefficient does not change by varying /3: 

a = 0.2142 ±0.0259 (6.21) 

This is again in agreement with previous analyses. On the other hand, if we admit a 
restoration of the U{1)a symmetry with temperature, we find a dependence on the 
(3 coefficient. 
For f3 = 0.2 we have 

a = 0.1484 ± 0.009 for Ti = 160 MeV (6.22) 
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and 



a = 0.1196 ±0.0161 for 



1 — 



100 MeV. 



(6.23) 



For (3 = 0.3 we have 



a = 0.1024 ±0.0042 for Ti 



1 — 



160 MeV 



(6.24) 



and 



a = 0.08101 ± 0.007 for Ti = 100 MeV. 



(6.25) 



However, according to Shuryak |134j it is very unlikely that the restoration of 
U{1)a can occur before chiral symmetry restoration; therefore, the value Ti = 
100 MeV should not be taken too seriously. In any case, it appears clear that 
restoration of U{1)a symmetry can strongly influence the behaviour of the criti- 
cal curve. In particular, the effect of the restoration of the axial symmetry with 
temperature could favour the splitting of the critical lines for the light quarks. It 
was pointed out in ^| that, for realistic values of isospin chemical potential, the 
strenght of the flavor-mixing interaction should prevent the separation of the critical 
lines. On this subject, the inclusion of the flavor asymmetry due to the different 
masses should be taken into account too. We are currently performing an analysis 
which considers Am = (m^ — m^) ~ m„ and a temperature-dependent coupling of 
the instanton induced interaction. 

Moreover, a study of the NJL model with an imaginary chemical potential is 
under investigation. This study would allow a direct comparison with the same 
analyses performed on the lattice. 

The mean-fleld NJL effective potential is symmetric for /i — ^ —fi, and due to its 
analytical properties depends on /i^ only, where here /i is a real chemical potential 
equal for all of the flavors. The critical curve depends on fi^ too. We can analytically 
prolungate /i ijJ'im to have that 



Therefore, in the regime of imaginary chemical potentials, we expect to flnd a critical 
curve which is the analytical prolongation to negative /i^ of our previous results. 



(6.26) 
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6.4 Conclusions 



In this application we have studied some general features of the QCD critical line 
in the framework of a NJL model. In section one, we have compared the physics at 
/iq ^ with that at /// 7^ 0. In section two, we have varied quark masses to show 
that the order of finite temperature transition changes if we consider small enough 
masses. In section three, we have studied the dependence of the slope of the critical 
curve Tc{ijl) on Nf. The fact that we have found a remarkable agreement between 
our analyses and lattice results may be taken as a proof that the NJL model can 
be trusted as a good toy model for describing the chiral transition. In addition, it 
is worth to stress that, since our analyses are performed at the mean-field level, we 
have basically computed the termodynamics corresponding to a gas of constituent 
quarks; on the other hand, Lattice-QCD analyses take into account bound states 
and fluctuations as well, therefore the successful comparison of the two approaches 
is even more impressive. 



117 



Conclusions 



In this work we have carried out a study of the phase diagram of QCD. We have 
started, in the Preface and in the first introductory Chapter, with a general survey 
on the QCD phase diagram. We have analyzed the phenomenon of chiral symmetry 
breaking and shown its connection with the confinement problem. We have discussed 
the high temperature regime, where dcconfinement should take place, and the low 
temperature and high densities regime, where color superconductivity is expected to 
occur. The presence of a critical point along the critical line has been highlighted. 

In the second Chapter we have introduced a subject which became rather pop- 
ular in the last few years and which motivated our work: the phenomenon of Bose- 
Einstein condensation of mesons in the low temperature regime of the theory. Ac- 
cording to several models, pion and kaon condensates are expected to form as the 
isospin and the hypercharge chemical potentials are high enough. This phenomenon 
is correlated to the "standard" QCD critical behaviour defined by the study of chi- 
ral symmetry breaking and restoration, because the non-perturbative QCD vacuum 
infiuences meson properties too. In this Chapter a review on this subject has been 
given. We have also presented some general aspects of QCD at finite chemical 
potential which are of interest for the present discussion. 

In order to carry out an original research project on this issue, we needed a 
model displaying the critical chiral behaviour and encompassing the possibility of 
describing meson condensates. Actually, due to the non-perturbative regime of 
QCD, these issues cannot be faced directly within the fundamental theory. A class 
of effective models reproduces a critical behaviour which is expected to be consistent 
with QCD. Amongst them, two microscopic models, ladder-QCD and Nambu-Jona- 
Lasinio, have been our tool to investigate the problem. In these models quarks are 
treated as fundamental degrees of freedom; chiral symmetry and meson condensation 
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are studied on the same footing, as a dynamical effect, by varying the thermody- 
namic parameters. We have hmited our analyses to the regime of moderate chemical 
potentials; superconducting effects have not been taken into account. 

In the third Chapter we have considered pion condensation in the framework 
of the ladder-QCD model. This model is able to highlight the critical behaviour 
related to the chiral symmetry, and gave a prediction on the existence of a critical 
point along the (T, fi) critical line from the late '80s on. In the present application, 
we have found that the critical value of isospin chemical potential fif for having 
pion condensation (at T = = 0) is just half of the pion mass. This result is in 
agreement with many other analyses obtained within effective models and lattice 
simulations. We have also studied the effect of an isospin chemical potential on 
the splitting of the critical lines relative to the light quarks. This fact could be of 
great interest in heavy ion collision experiments, in the attempt of finding a definite 
signature for the transition to the deconfined phase of the theory. 

In Chapters four and five, we have studied the subject of pion and kaon con- 
densation by using the NJL model. The use of a model different from ladder-QCD 
was motivated by the technical difficulties encountered for extending the model at 
finite (T, //), when considering meson condensates. The NJL model is simpler to deal 
with, and encompasses the relevant features for describing the problem as well. The 
results we have found confirm the possibility of having pion and kaon condensates 
in the regimes of high isospin and hypercharge densities, respectively. In particular, 
in Chapter four we have performed an analysis of the two flavor model in the whole 
space of thermodynamic parameters {T, Hg, Hj). In Chapter five we have extended 
the study to the third flavor, and studied the competition between pion and kaon 
condensates in relation with the behaviour of usual chiral condensates. We have 
proposed an original hypothesis for kaon condensation, by considering a chemical 
potential associated with light quarks only. We remark that we have neglected in 
these works any superconducting effect; a natural extension of this analysis would 
require the combined presence of di-quark and meson condensates. 

In Chapter six, we have focused our attention to the critical line for chiral symme- 
try restoration, by neglecting meson condensates. We have studied the dependence 
of the critical temperature on quark and isospin chemical potentials. We have 
shown that the order of the zero chemical potential transition changes from second to 
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first order if we consider low enough masses and a coupling between fiavors induced 
by instantons ('t Hooft). We have studied the dependence of the critical curve on a 
chemical potential, by varying the number of massless fiavors Nf. This study within 
a microscopic model, here the NJL model, is motivated by the growing interest that 
the lattice community has shown for this issue; actually, in recent years many pos- 
sible algorithms have been proposed to overcome the difficulty (the sign problem) 
of the fermion determinant at finite baryon chemical potential. Our results are di- 
rectly comparable with Montc-Carlo simulations. A study of a Ginzburg-Landau 
expansion of the free-energy, in order to highlight the cubic term induced by instan- 
tons, is under examination. Another development regards the study of the critical 
line as a function of an imaginary chemical potential. The effect of axial symmetry 
restoration and of the explicit isospin breaking (Am = — m„) on the behaviour 
of the light quarks critical lines is under study as well. 
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